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Abstract 
records is treated for the case when the input and output quantities have random perturbations. 


The problem of determining process dynamic characteristics from process operating 


The correlation and spectral analysis techniques necessary are briefly reviewed and then applied 
to some data obtained from a computer simulation of a stirred tank reactor. 
was used for analysis 


\ digital computer 
; and the effect of length of record, the spacing between points, and type 
of weighting function used to estimate the process transfer function and impulse response are 


considered. The effect of process non-linearity is also considered. 


Résumé— On traite dans ce travail le probléme de la détermination des caractéristiques dynami- 
ques d'un systéme a partir d’enregistrements en cours de fonctionnement dans le cas of “ signal ” 
et * réponse ** sont soumis a des perturbations aléatoires. Les techniques d’analyse de la corréla- 
tion et d'analyse spectrale sont rapidement passées en revue et appliquées a des données tirées 
d'un caleulateur simulant un réacteur a agitation. On a utilisé pour lanalyse un calculateur 
digital (ou numérique). On étudie leffet de le longueur de lenregistrement, de la distance entre 
les points, et le type de fonction “ poids,” utilisé pour estimer la fonction de transfert et la réponse 
impulsionnelle. On envisage aussi l'effet de la non linéarité du systéme. 


Zusammenfassung—Die Ermittliing der dynamischen Eigenschaften eines Prozesses aus der 
Registrierung der Arbeitsweise des Prozesses wird hier fiir den Fall behandelt, wenn Einsatz 
und Produkt zufalligen Mengenschwankungen unterliegen. Die Korrelation und die notwendige 
spektrale Analysentechnik werden kurz mitgeteilt und dann auf einige Werte aus dem Simulator 
eines geriihrten Tankreaktors angewendet. 
benutzt. 


Ein Digitalrechner wurde fiir die Untersuchung 
Der Einfluss der Registrierlinge, des Abstandes zwischen den Punkten und der Art 
der Gewichtsfunktion zur Bestimmung der Cbertragungsfunktion des Prozesses werden betrachtet. 
Auch der Einfluss der Nichtlinearitét wird beriicksichtigt. 


INTRODUCTION output are used here synonymously with cause 


A KNOWLEDGE of process dynamics is important 


in the design and analysis of processes and their 
control systems. 
dynamics of a process might be obtained from a 
detailed analysis of the basic phenomena or they 
might be obtained from observation of an operat- 
ing process or model of the process. The latter 


The equations describing the 


method will be treated in this paper. In brief, 
the dynamic properties are deduced from simul- 
taneous observations of the inputs to the process 


and the outputs from the process. Input and 


and effect. Typical inputs to a chemical process 
might be reactant concentration, flow rate. or 
valve position. Typical outputs might be product 
concentration, temperature, or viscosity. The 
treatment presented here assumes that one or 
more uncoupled, randomly varying inputs are 
present. The desired dynamic property is the 
transfer function or impulse response relating 
one of the inputs to the output being considered. 
A digital computer was used for the computa- 
tions required, 
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In the next section the basic relations for 
analysis of linear systems are briefly reviewed ; 
following sections some data obtained 


stirred tank 


in the 
from a mathematical model of a 
reactor is used to demonstrate the application 
of these methods to a non-linear process and to 
determine good values for some of the parameters 


used in the analysis, 


THEORY 


If the input toa Process 1s 4 (?) and the output 
is y(t), then for a linear process the output ts 


given by 


y (t) h,,(r)a(t r)dr (1) 


-@ 
where h,, (7) is the impulse response or weighting 
function for the process and is one means of 
expressing the dynamic characteristics of the 
process. Alternatively the dynamic properties 
can be expressed in terms of frequency as follows : 

x 
Y (w) 


” 
X (w) (2) 


H (w) h (+r) exp (— jwr) dr 
-@ 

where H (a) is the frequency response or transfer 
function of the process and is the Fourier trans- 
form of h,, (r). Y (w) and X (w) are the Fourier 
transforms of y (t) and 2 (t), assuming these trans- 
forms exist. Ifa (t) is an impulse or step function, 
h,,, (+) can be determined directly from equation 
(1). If a (t) is a sinusoid, H (w) can be determined 
most easily from equation (2). In general, equa- 
tions (1) and (2) can be used to determine the 
dynamic properties if x (t) is the sole cause of the 
variations in y (t) and each is known without error. 
The latter restriction is not too limiting, but the 
requirement that there be no other sources of 
output variation than the input being studied is 
a serious practical limitation to direct Fourier 
analysis of the time series. To remove this 
restriction the concept of correlation between 
input and output signals must be introduced. 
The correlation function is defined as: 


T 
q 


-T 


dy (7) = lim 
T > 


a 
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where ¢ is the autocorrelation function if f and g 
refer to the same signal and the cross-correlation 
function if they different From 
equations (1) and (3) it can be shown that 


are signals, 


h,, (72) bax (tT — 73) dry (4) 


The Fourier transform of the correlation function 
[®,, (w)] is called the spectral density, and it can 
be shown that 


®,, (w) = H,, (w) On (aw). (5) 
It then follows that 
H,, (w) ’ (6) 


which is identical with equation (2) if the auto- 
correlation function of the input is considered 
the input to the process and the cross-correlation 
function as the output of the process. 

The primary difliculty in using equation (6) 
is in obtaining accurate estimates of the spectral 
densities. If a digital computer is used for 
analysis, a series of numbers must be used to 
represent the original time series. Since the only 
practical methods of analysis require equally 
spaced points in time, the # and y signals must be 
represented by a finite number of equally spaced 
points. 

Assume a(t) is available only for the times 
t=<k, At; k,=0,1,2...N, and y(t) for 
t = k, At, where k, M; M and N 
are finite; and the time reference is identical 
for w(t) and y(t). M need not equal N, but it is 
a practical choice and will be assumed here. 
The correlation function is then given by 
yg (7) ~ bjg (vAt) 

N—|s 
N 4 \v| ya Pe lh To , 
+ |r| At); OY <N 


1 N 


: y k, At + 
N n | = IY] ay 7 T 


+ |v] At) g(k, At); —N 


where an asterisk is used to indicate an estimate. 
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Similarly equations (4), (5) and (6) are replaced 
by 

N 

dy (vAt) 2 


¥g=0 


?, (AAw) 


hy, (vg At) 2, (vAt — vg At) At, (8) 


H?, (Adw) ®, (AAw) (9) 
and 
?, (AAw) 


4 . (10) 
®, (AAw) 


H*, (AAw) 
The spectral densities are obtained from the 
correlation functions as follows : 
N 


@7, (AAw) 2 W (vAt) dy, (vt) 
N 


[cos (AAwwAt) — j sin (AXowAt)] At (11) 
where W(»At) is a weighting function and will 
be discussed later. 

Finally the impulse response is obtained from 
H}, (Adw) by 


he, (vAt) Hy, (AAw) 


[cos (AAwwAt) + J sin (AAwv Mt) | (12) 

Since data is available only at discrete intervals, 
any frequencies higher than 1,2 At cycles sec 
will not be detectable. However, this information 
is not ignored but is folded back to the lower 
frequencies giving an erroneous estimate of the 
power in the lower frequencies. This phenomenon 
is called * aliasing * and will be discussed later. 
The fact that finite amounts of data are used 
implies that more data is used to compute 
yj, (v At) for small |»| than large |v}. 
reliable values must thus be weighted lower than 


The less 


the more reliable values. 

If other random disturbances are superimposed 
on y (t) it can be shown that these same relations 
are valid provided the disturbances are not cor- 
related with a(t). A measure of the amount of 
extraneous signal present is given by the coherency 
of the signals defined as 


Dry (w)? ; 
PD, (w) o,, (w). 


A coherency of 1-0 would indicate no noise 
present. A coherency of 0 would correspond 
to no correlation at all between the signals. 


1 


More detailed of the basic 
principles involved in spectral analysis can be 


found in references [3], [4], [5] and [6]. 


presentations 


Process SIMULATION 
To demonstrate the use of spectral techniques, 
to investigate the effect of non-linearities on the 
of to 
determine suitable values for analysis parameters, 


estimation process characteristics, and 


a mathematical model of a stirred tank reactor 


was used to generate data. The process model 


has been described in detail previously [1], [2], [7 

and is only summarized here. 
The 24 

Heat takes place in a stirred tank reactor. Heat 


irreversible reaction > Products 


is removed from the reactor by a cooling coil. 

The detining equations are : 

dA, 
dé 

_ aT, 

. dé 
kV A®, (AM) 


2He Pe Ce 
hA,, 


J 


@ (A, Ay) 


| p Qp C ( T, T,) 


hA,F 


ky exp (16) 


E 
| RT, 


If these equations are linearized about the 
the 
differential equation relating outlet concentration 


steady state conditions, following linear 


to init concentration and inlet temperature can 
be derived : 


da day 


2 

ee tee 
ay aC" di 
The 


in another paper [7], though some of the numerical 


Byag ~ Bsty (17) 


expressions for the constants may be found 


values were changed for this study. 

A model of the system was programmed for 
the Univae Scientific 1103 Computer. Inlet 
temperature alone and then inlet temperature 
and inlet The 
output used was The 
inlet disturbances were created by generating 
random numbers with a 
standard deviation of unity. 


compositions were disturbed. 


outlet concentration. 
mean of zero and a 
These random 
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numbers were then multiplied by a factor to 
create the desired size of disturbance and added 
to the steady-state value of the input to obtain 
the input actually presented to the process. 

The simulation programme used the Runge- 
Kutta Gill Method of integration employing a 
time interval of 8 sec. The random disturbances 
were introduced as step changes occurring at 
intervals of 32 sec, 128 sec, or 256 sec. Ease of 


manipulation and computation restricted to 


H, ,(w) [jw (1 32 10 3) (6 S4 


and for run 1 


i, (w) [jw (3 a3 


The latter run was included to demonstrate the effects of changing process parameters, 


10 *)| [jw 


W. Treeney 


901 the number of points which were handled as 
one continuous piece of data. If longer runs were 
desired, values of input and output were recorded 
in several sets of 901 points each. Thus a run of 
9010. points was stored in ten sets of 901 points 
‘ach. Run conditions are summarized in Table 1. 
Input and output quantities were stored for each 
time point. 
For runs 2 through 5 the linearized transfer 


function was 


06-0501 
10 J] [jw 


(1-32 . 10%) — (6.84 « 10%))] 


0-636 
(4-56 10 *)] 


In the 


former process the controlling time constant is 759 sec and in the latter 219 sec. 


Table 1. 


Run conditions for process simulation 





Run T 


(see) 


- 
72, 


144,160 


144,160 


Oso 


Output 
concentration 


Input Input 


step duration temperature 


(sec) variance variance 


5-08 1-* 


1-03 104% 
103 « 10°34 
wr? 


6-20 #38 « 10 


15-50 506 © 10 


16-31 102 « OF 





*Run 1 was made using the short time constant process, all other runs used the long time constant process, 


tRun 2A was obtained by rearranging the data of Run : 


tRun 5 was made using two input disturbances. The second disturbance was a stepwise variation in inlet concentration 


with random amplitude and 128 see duration, Variance was 3-0 


COMPUTATION OF TRANSFER FUNCTION 


Transfer functions were calculated as follows : 
1. The mean value of 
determined by using all values recorded 
for that variable. 
the generated input values, 
This mean was then subtracted from all 
the individual values of the variable to give 


each variable was 


This was done even for 


a time series with a zero mean. 
The then correlated by 


* sets” 


variables were 


using equation (7) and the “ set 


1-3, 


correlations " averaged to determine the 
“average correlation.” Some 
information is being discarded by ignoring 
the products of points at the 
junctions of consecutive sets. 

The 
weighted. 

Spectral density functions were computed 
from the average correlation 
functions using equation (11). 


potential 
lagged 
functions were 


average correlation 


weighted 


The transfer function was determined from 
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Table 2. 


disturbances 


Parameters used for spectral analysis 





T» 


(sec) 


1,024 


2.048 
8.192 
16.384 


2.048 


8,192 
4.006 
16584 


£006 
6.144 
8,192 


Mt T 


€ 


0-0366 


O-O105 
O-O422 
O-O422 


00-0422 
0-422 


O-O422 


0-0422 
0-422 
O-O422 





the 
the spectral density of the input. 


The impulse 
determined from the transfer function. 


ratio of the cross-spectral density to 


response if desired was 


Steps 4 and 5 can be reversed by computing the 
spectral density of the correlation functions and 
convoluting this with the Fourier transform of 
the weighting function. The methods are equiva- 
lent and 
convenience. 

A number of different weighting functions can 
In this work three The 
first, which is the simplest possible, is the trunca- 


the choice depends on computing 


be used. were used. 


tion estimate, defined as follows : 


Lb (vt) 


exp (— jAAwvAt) At; m<N 


1®;, (AAw) 


(18) 


The estimate is Bartlett's estimate: 


second 


pPy, (ALw) 


E (1 fl) steam 


exp (— jAAwvAl) At; m< N (19) 


The third weighting function (called ““Hamming”) 
is given by 
nz, (Aw) 
+ 0-547 BF, (AAw) + 0-237 DF, [(A 


0-237 , [(A — 1) Aw] + 
1) Aw] (20) 


Using the first method the lagged products, 
although available to NAt, are treated as zero 
beyond mAt. In the second method the lagged 
products are also treated as zero beyond mAt 
but are linearly reduced between 0 and mAt. 
The third is essentially a procedure for smoothing 
the truncation estimate. In general the truncation 
estimate has less bias and more variance than the 
Bartlett estimate [5]. 

The primary computational variables are thus 
the sampling interval, At; the largest lagged 
mMt=T,,; total 

and the number 


product, run 


NAL = T,; 
obtain the average correlation function, P. 
Table 2 
computation are summarized. 

The sampling affects the 
frequency recorded and hence the ~ aliasing.” 
The total run length affects the stability and 


length, 
used to 
In 


the values of the parameters used for 


of sets 


interval highest 


accuracy of the estimates, a twofold improvement 
fourfold 
increase in total run length. Thirdly, the highest 
lag used will determine the frequency resolution 
which can be achieved. 


in stability or accuracy requires a 


Discussion or RESULTS 

The results of the computations are presented 
in Figs. 1 through 7. In each figure a solid curve 
is used to present the results obtained using 
the linearized equations. Deviations are due to, 
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first, the errors of estimate caused by the limited 
amount of data available and, second, the fact 
that the data 


pre CSS, 


was obtained from a non-linear 


Input spectral density 


In the frequency range of interest, it is impor- 
tant that the input spectrum be relatively flat 
or white. Obviously, if no power were present 
at some frequency in the input, there could be 
no information concerning the process behaviour 
at the frequency. Serious analysis problems can 
Because of lack of 
resolution in the input and output there will be 


develop when this happens. 


some “ smearing ” of the spectrum whenever the 


spectrum has large variations and erroneous 


estimates will be obtained in such a region. 


e 
» 
> 
e 
> 
3 
x 
é 
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closely followed this theoretical spectrum and are 
not shown. 


Variation of T,,, T, and T,, 


The accuracy of the estimates obtained from 
an analysis is determined by the parameters, 
a] 

T..» T.. T, and At. 
these parameters it is necessary to have some idea 
of the frequency spectrum of interest in the 


Before selecting values for 


process. This can be expressed by designating an 
upper frequency and a frequency interval. The 
latter the resolution fre- 
quencies or the low frequency response of the 


determines between 
process, One means of characterizing the frequency 
spectrum of interest in a process is by means of the 
the largest significant time constant, here called 
T.. Many processes will not pass significant 


SE 
00 








wo al 


wh | 


Theoretical spectral density of input showing relative amounts of power at important frequencies 
when duration of random step is four sampling intervals. 


When this 


whitening may be helpful. 


occurs, techniques such as pre- 
It is also important 
that the amount of power in the spectrum be 
In Fig. 1 
is shown the theoretical spectrum for one of the 
temperature inputs used in this study. This 


cos 


small beyond the folding frequency. 


spectrum, which has the equation Ko*p (1 
yw) ‘w*, is smooth in the range of interest, 0 to 
and contains little power beyond the folding 
frequency, ,. However, there are some points of 
little or no power, which caused some difficulties 
in taking inverse transforms as will be discussed 
later. The spectra as determined from the data 


Was 


energy at frequencies more than two decades 
above 1,7’, rad /sec, and resolution to 1/7, rad /sec 
is often sufficient. By estimating 7, an approxi- 
mate frequency range of interest can be estab- 
lished. 

The analysis parameter that determines the 
highest frequency about which information will 
be available is At, the sampling interval. 1, 
determines the resolution possible. 7’, and 7’, 
must be selected sufficiently large that statistically 
valid information is obtained. These points are 
demonstrated by considering some results obtained 
using different combinations of parameters. 
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In Table 2 it will be noted that two processes 
were studied— one having a T, 3-47 times the 
other. In each case the At used was such that the 
highest frequency studied was well above the 
point at which the process transmission had 
declined to a negligible amount. The upper 
frequency limit is thus adequate for all cases 
studied. 

The effect of varying T,, and consequently the 


am 


resolution is shown in Fig. 2, which presents the 
results of runs 2, 2Aa, and 2Ab. In run 2, 7,,/7, 
is too small (2-7) to give sufficient resolution, 
and severe damping results. Run 2Aa, which is 
almost identical except for 7',,/T, (10-8) is a 
much better approximation to the linearized 
transfer function. This same point is also shown 
by comparing run 2, Fig. 2, and run 1, Fig. 4. 
Although a longer T,, 


_ is used in run 2, the results 


IMAGINARY PART 





W* RADIANS/ W 
ws 333.3 SECONDS 








REAL PART 


| Ty * 2048 SECONDS, RUN 2 


2 "y = 8192 SECONDS,RUN 2Aa 


3 T,, = 16384 SECONDS RUN 2Ab 


4 LINEARIZED 


Transfer function of the long time constant process showing the effect of varying T,, relative 
to T,, Tp and T,. 
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IMAGINARY PART 








REAL PART 


1. LINEARIZEO 
2 HAMMING 


3 BARTLETT 
© TRUNCATION 





Transfer function of the long time constant process. The difference between types of estimates is shown 


Fie, 3. 
160 





Determination of dynamic characteristics of processes in the presence of random disturbances 


are not nearly as good as in run 1 because 7’, has 
decreased from 759 sec in 2 to 
run 1. In general, both runs 1 and 2A gave 
reasonable approximations to the linearized 
transfer functions. Part of the difference, it must 
be remembered, is due to non-linearities in the 


run 219 sec in 


process. Run 1 has a much smaller disturbance, 
which accounts for part of the improvement in 
the estimate. 

In addition to selecting a proper Af and T7,, it is 


necessary to have enough data to obtain good 


H — 
1a 


estimates of the parameters. This implies that 
T,, and T,, must be of proper length. Comparison 
of run 2Aa and run 2Ab, Fig. 2, shows that 
increasing T,,, to avoid damping is counteracted 
by loss of stability if 7, becomes too large 
relative to T,, and T,,. The results of runs 4a, 
th, 4c, 5a, Sb and 5e, 
illustrate this point. 


while not shown, also 
Runs 4a and 4b both give 
satisfactory results, while 4c, in which T,, is too 
large relative to T,, and T,,, exhibits poor stability. 
Runs 5a, 5b and 5c all give satisfactory results, 


IMAGINARY PART 





| 


10.0 W=RADIANS/Y (333.3 SEC) 














Fic. 4. Transfer function of the short time constant process. 


| LINEARIZED 
2 EMPIRICAL 


Parameters unsatisfactory in analysing the 


long time constant process are satisfactory here. 
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showing that 7, can be varied by a factor of two 
and still give satisfactory results. Thus, in order 
to be able to use a high value of T,,,, a correspon- 
dingly large T,, and T,, must be used. In general, 
with a reasonable band width process, one to two 
thousand points are necessary as an absolute 
minimum, Often more will be 


many points 


desired. The only practical way of performing the 


computation is by means of a digital computer 
such as the Univac Scientific Model 1108 used in 
this work. 


Type of weighting function 
Fig. 3 presents the results obtained using three 
different weighting functions on the data from 
run 4. No curve is drawn through the points for 
the the 
fluctuations. In general the Bartlett and Hamm- 


truncation estimate because of large 


ing estimates are smooth curves and good approxi- 
mations to the transfer function. These points 
are not shown but are at the same frequencies 


as the truncation points. The Hamming has 


Hlowan and J, W. Trerney 


somewhat less bias but a higher variance than the 
Bartlett. 
be used, but these two are particularly simple to 
the Bartlett in the time domain and the 


Many other weighting functions could 


apply 
Hamming in the frequency domain. For all the 
other results reported here the Bartlett weighting 
was used, primarily because the heavier smoothing 
gave better impulse response curves as discussed 


below. 


The impulse response 
The impulse response is obtained from the 
transfer function using equation (12), and is an 
alternate way of describing the process dynamics. 
The from 7 Atrad sec to 
m7 Atrad sec was used for the inverse trans- 
procedure for 


entire range 


form without any weighting. 


Since the input disturbance contained little or no 
power at some frequencies covered (Fig. 1), the 
transfer function was unreliable at these fre- 
quencies. In those cases where the errors resulted 


in estimates of high transmission at these high 


2900 











Fig. 5. 


2048 SECONDS, FUN 2 


Tw * 6192 SECONDS, FUN 240 
Ty * 16384 SECONDS, RUN 2AD 
LINEARIZED 


™w* 


Impulse response of the long time constant process. The damping effect of 7,, is shown here also. 
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frequencies, taking the inverse transform without 
the 
estimates of 


resulted in 
the 
Bartlett estimates of the transfer function were 


smoothing values 


incorrect 


obviously 
impulse response. 
less susceptible to this kind of error and were 
used to compute the impulse response. 


Fig. 5 shows the impulse response curves as 


estimated from the transfer functions of runs 2, 
2Aa, and 2Ab, which were shown in Fig. 2. The 
effect of the choice of T,, is as obvious here as 
for the Fig. 6 
effect of attempting to weight the results of Fig. 5 


transfer function. shows the 
by the factor m/(m 
Bartlett weighting. 


but it is definitely not an adequate remedy for a 


v), Le. the reciprocal of the 
Some improvement results, 


poor choice of T,. 


The Effect of Non-linearities 
Fig. 7 presents the results obtained from runs 
made with disturbance. 


varying severity of 


The variance of the output concentration which 


TIME ~ SECONDS 


resulted from these disturbances is shown as the 
parameter for these curves. It may be noted that 
as the variance increases, the impulse response is 
This 
agreement with the results of other investigations 
Run 


3 differed from run 2 only in the degree to which 


more damped. greater damping is in 


of non-linear processes of this type [1, 7]. 


the process was disturbed and the only detectable 
difference in the results existed in the peak height 
as shown in Fig. 7. Runs 4a, and 4b gave essen- 
tially identical results as did runs 5a, 5b and 5e ; 
consequently only single curves are shown for 


runs 4 and 5 in Fig. 7. These results indicate 


that satisfactory results can be obtained using 
a wide range of parameters in the analysis. Run 
te which used too low a value of 7/7, exhibited 
It should be 
noted that the largest single disturbance used, 


poor stability and is not shown. 


run 4, produced an outlet concentration which 
0-071 moles /ft*. 
The magnitude of this variation can be better 


had a standard deviation of 








Fig. 6. 


Adjusted impulse response of the long time constant process. 


A factor m/(m — v) applied to the 


estimates of Fig. 5 improves the results but can not compensate for the effects of poor analysis parameters. 





TIME — SECONDS 


L 1 
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Fig. 7. 


Impulse response of the long time constant process. 


The non-linearities of the process show a damping effect. 


appreciated by considering that the concentration 
was outside the range 0-583-0-817 moles ft® more 
than 10 per cent of the time. Since the desired 
value is 0-700 moles ft® from a physically possible 
range of 0-1-0 moles ft® such fluctuations would 
be intolerable in any operational unit and are 
well in the non-linear range. In spite of this the 
linear approximation, given by the solid curve, 
and the calculated curve for run 4 are in good 
agreement. To the this 
can be considered representative of chemical 


extent that process 
processes then it can be said that approximations 
as obtained by this method are satisfactory for 
most engineering purposes. 


6. Multiple Inputs 


Inruns 1 through 4 only one process disturbance 
(inlet temperature) was present. 
to make the interpretation of results simpler. 
In practice this would rarely be the case, and, in 
fact, for a single input and no random measure- 
ment errors a conventional Fourier analysis 
would be sufficient. In run 5, Fig. 7, is shown the 


This was done 


impulse response computed for simultaneous 


and inlet 


uncorre- 


disturbances in’ inlet temperature 


concentration. The disturbances were 
lated, although their spectra were similar. As 
indicated previously, the techniques described 
in this paper are valid for multiple inputs only 
if the inputs are uncorrelated, The techniques 
can be extended to correllated multiple inputs 
by solving simultaneous equations but this is 
beyond the scope of this paper. The additional 
input disturbances resulted in larger excursions 
of output concentration (a, 0-1) than for 
run 4. In spite of this a satisfactory approxima- 
tion to the impulse response was obtained. The 
coherency for run 5 was calculated to be about 
0-43. 
for run 4 in the frequency range where significant 


By comparison, the calculated coherency 


transmission occurs was about 0.96. 


CONCLUSIONS 


The techniques described in this paper have 


been shown to give good approximations to the 


dynamic characteristics of a typical chemical 
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process when random uncorrelated disturbances 
to the process are present. Although the method 
uses an essentially linear technique to study a 
non-linear process, the results seem adequate for 
most engineering purposes. A digital computer 
can be efliciently used to perform the calculation. 

In using the technique care should be exercised 
at all stages of the analysis to insure that no 
serious aliasing occurs and that suflicient resolu- 
tion and stability are achieved. “ Aliasing” of 
the input should be avoided and the input should 
If an 
approximate value of the largest significant time 


be nearly white in the area of interest. 


constant of the process, T’,, is available, the results 
of this study indicate that a good analysis should 
be achieved by using At ~ 0.04 T., T,, ~ 10 T., 
and T,, > 100 T.. If the data is to be gathered 
- 35 T, and it would 


Approxi- 


in pieces, T, should be 


, 
be desirable to increase T,, slightly. 
mately 2500 points are required by these values 
which should provide sufficient accuracy and 
resolution in the important frequency range. 
If * aliasing ” 
input, sampling the input at a higher frequency 


Is suspected in the spectrum of the 


and comparing spectra will help in determining 
if this effect is serious. If the power of the 
input is not constant, or nearly so, at the lower 
frequencies, pre-whitening the spectra may be 
advantageous. 
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NOTATION 
dimensionless concentration 
heat transfer area ft? 
steady state value of concentration Ib mole/ft® 


concentration of outlet and inlet stream 
Ib mole ft® 


hr) 
II (@) 
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heat capacity of process stream, coolant 
B.t.u/lb °R 
B.t.u/lb mole 
B.t.u/sec ft? °R 


stream, and tank contents 
activation energy 

heat transfer coeflicient 
impulse response of process 
transfer function or frequency response of 
process 

ft? /lb mole sec 
ft? /Ib mole sec 


reaction velocity 
reaction velocity constant 
number of lag intervals used in Truncation 
or Bartlett Estimates 
number of time intervals used in output 
and input data 
number of time intervals represented in 
each piece of data 
number of pieces of data 
inlet and outlet flow rates ft® /hr 
B.t.u/ Ib mole °R 
time, sec ; dimensionless teraperature 

TT, 
largest time constant of process 
inlet coolant temperature 
highest lag used in taking transforms 

mM 
total data length — NAt = PT, 
length of each piece of data pat 
steady state value of temperature 
temperature of outlet and inlet streams R 
volume of tank ft? 
input and output of a process 
length of step used in disturbance sec 
heat of reaction B.t.u, lb mole 
sampling interval sec 
frequency interval rad unit time 
time sec 
integer used with frequency interval to denote 


yas constant 


frequency 

integer used with sampling interval to denote 
lag 

Ib /ft8 
Ib /ft® 


density of reactant mixture 
density of coolant 

standard deviation of a variable 
time lag 

correlation function 

spectral function 


= dimensionless time qWy 


frequency rad unit time 
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Heat transfer studies between particles and liquid medium in a 
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Abstract—The heat transfer coeflicients between the particles and the fluidizing medium under 
unsteady state conditions, are determined by transient heating of cold water stream, with hot 
solid particles dropped into the fluidizing bed, in a specially designed liquid fluidizing column. 
Heat transfer coeflicients varying from 113 to 620 B.t.u hr ft®) F are obtained and these values 
increase with increase in particle diameter and mass velocity. The present data are also compared 
with previous published data on gas fluidized beds. 


Résume—Les coeflicients de transfert de chaleur entre les particules et le milieu fuidisant 
dans des conditions d'état non stationnaire sont déterminés par le chauffage fugace du courant 
d'eau froide, & Taide de particules solides chaudes projettées dans le lit fluidisé dune colonne 
de fluidisation liquide. 

Les auteurs obtiennent des coeflicients de transfert de 113 a 620 
B.t.u/hr ft®/°F et ces valeurs augmentent avec laccroissement du diamétre des particules et 


chaleur variant de 
la vitesse massique, Ces données sont comparées avec celles publiées antéricurement sur les 
lits de gas fluidisés. 


Zusammenfassung—Die Wiirmetibergangskoeflizienten zwischen den ‘Teilchen und dem 
fluidisierenden Medium im instationiren Zustand wurden dadurch bestimmt, dass cin Valter 
Wasserstrom durch in das Fliessbett geworfene heisse Teilchen in einer dafiir besonders cingeri- 
chteten Kolonne nichtstationir erwiirmt wurde. Die Wirmeitibergangskoflizienten betrugen 
$40 bis 3500 W/m?® grd und stiegen an mit dem Teilchendurchmesser und mit der Massenges- 


chwindigkeit. Die erhaltenen Werte wurden mit friher verOffentlichten Daten fiir Gas-Fliessbetten 


verglichen. 


Struptes on heat transfer between solid particles 
and a fluidizing medium are necessary to clucidate 
the mechanism of dissipation of heat of reaction, 
heat of dilution etc., in a fluidized bed and also 
in fluidized bed reactor design. Kererrenrine 
et al. {1} found the heat transfer coeflicients 
between the particles and the fluidizing medium 
to vary between 3 and 9-1 B.t.u/hr/ft.2 °F and 


explained that the low values of the coetlicients 


were due to the reduction in the slip velocity of 
the gas past the particles. An attempt was made 
by Watton et al. [2] to determine the temperature 
But it 
was not clear whether the thermocouple indicated 


gradient between the gas and the solids. 


the true temperature of the solid, or of the fluid 
or of some mean between them. Ercunorn and 
Wuirr [3] were unable to measure the temperature 


difference between the solid and the gas although 


16 


the particle temperature was maintained constant 
by dielectric heating. 

The 
avoided by Wams.iey and Jonuanson [4] by using 
How- 


estimation of solid temperature was 
unsteady state experimental conditions. 
heat transfer coetlicients were 
0-06 0-9 


short 


ever, very low 
obtained, ranging between and 


B.t.u,hr ft?,°F. This the 
circuiting of the gas, in the form of bubbles and 


was due to 
channels, without attaining thermal equilibrium 
with the bed. Recalculating the 
Wams.ey and JoHanson, Suen and JounsTone 
[5] found that the heat transfer coefficients in 
discontinuous phase increased with gas velocity 


results of 


contrary to Wamsiey and Jonanson’s observa- 
tion. 

In general it is observed that very low heat 
transfer coeflicients were obtained in gas fluidized 
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beds [1, 2 and 4], which were mainly due to the 
formation of gas bubbles and short circuiting of 
gas through channels. However, in liquid fluidized 
the bubble channelling 
phenomena are absent. heat 
transfer coeflicients can be obtained if the solid 


beds formation and 


Hence, higher 
bed is fluidized by a liquid medium. To date, no 
data on such a system have been reported. 

The theory proposed by Wams.ey and Jouan- 
son, with minor modifications, has been used in 
this The 
because the direction of heat transfer was from 
the solids to the fluid, instead of from the fluid 


work, modifications were necessary 


to the solids. From the temperature-time history 
of the outlet liquid, heat transfer coeflicients 
were calculated from the following equations ; 


(t t;) a exp (b@) 


where a=(t 


h, A uw C, 
WOs(h, A 


and b 


° 
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and liquid medium in a fluidized bed 


EXPERIMENTAL APPARATUS 
PROCEDURE 


AND 


The apparatus used is shown in Fig. 1. Water was 
pumped at a constant flow rate and at a constant tempera- 
ture through the 3-25 in. LD. fluidizing column B, and 
also through 2-675 in. I.D. pressure equalizing column, C, 
where a known amount of solid particles (1000 ¢) were 
fluidized. Another batch of closely-sized particles equal 
to the particles present in the column, C, were carefully 
heated to a uniform temperature and after ascertaining 
whether the particles were at a constant uniform tempera- 
ture, they were transferred to the solid dropping device 
(Fig. 2) and the vessel was immediately closed by tying 
ona filter paper, N. The device was fixed to the apparatus 
and X”, the gate valve, K, being kept 
closed, Then the gate valve, K, quickly opened and solids 


at positions X 


were dropped into the fluidizing column, B, by puncturing 
the filter paper N with the steel rod M, and immediately 
closing the glass expander top by rubber cork Y. The 
outlet water temperature, read by a long-stem Beckmann’s 
thermometer f), having a range of 5-5 C, was noted as at 
function of time. This was done by starting a series of 
eight stop-watches simultaneously and stopping each at 
marked regular steps of temperature rise in the outlet 
water. 



































Fic, 1. 
A Solid holding chamber 
D Rotameter 
b, «, K Gate valve 


Beckmann thermometer (at fluid Outlet) 


a, d, e, f, wg. hy, Valves 


to 


B Main fluidizing column Cc 
S Water storage tank 


Experimental apparatus. 


Pressure equalizing column 
P Pump 
t, thermometer (at fluid Inlet) 
(Pipe-lines used are 1), in. standard size) 
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Fic. 2. Solids dropping device. 
Solid holding chamber PP’ 
Gate valve 

Steel rod x’, 
Filter paper 

Cylindrical vessel Y. 


Cork assembly 

Glass expander 

X” Fixing points of 
assembly 

Y’ Rubber corks 

Z Copper tube 


The variables studied were as follows : 

1. Particle sizes of : 10-14, 14-16, 16-20, 20-24, 24-28 
and 28-32 standard Tyler mesh, for granite and 
quartz materials. 

2. Mass velocity of water from 12,000 to 20,000 Ib/hr/ft®. 


2-205 Ib (1000 g) of solid material for exch run. 
The bed height in the main column, B, was veried 
between O-7ft and 1-5ft. 

The inlet temperature of the water was varied 
between 70 F and 85 F. 

The range of solids temperature, before transferring 
to the solids dropping device was varied between 
195 °F and 210 F. 

The specific gravities were determined as : 
2-66; quartz, 2-59. 

The specific heat of granite was found to be 0-21, 
and of quartz 0-18. 

The shape factor for the irregular solids was fixed 
at 1-75 [6]. 


uranite, 


Tyrpicat OBSERVATION 


DISCUSSION OF 


AND THE 


rue Resvuuits 
From the noted outlet temperature of the water 
from the thuidizing bed, the temperature differ- 
ences between the outlet and the inlet water of 
the fluidized bed (t 
time, #, on semi-log paper. It is observed that the 


t;), were plotted against 


middle portion of the temperature-time history 
curve resulted in a straight line on the semi-log 
paper. This 
with Wams.ey and Jonanson [4]. In this period, 


observation was in accordance 
the accuracy of the experiment is quite high; 


during the initial period of experiment, the 


disturbances in the flow pattern due to the 


dropping of the solids result in the erroneous 
observation of the outlet water temperature and 
in the last phase of the experiment, the mixing 
of hold-up water from chamber A (Figs. 1 and 2) 
above the outlet of the main column with the 
incoming water, results in an incorrect outlet 
water temperature. 

In general, for each set of experimental condi- 
tions, the same experiment was repeated at least 
three times to obtain reproducible results. In 
many cases reliable and reproducible values of 
heat transfer coetlicient were obtained. A typical 
observation of such a case is given in Table 1 for 
run Q.W.L.1. The plot of In (¢ — t,;) vs @ is shown 
in Fig. 3. The 6 values were found to be 2060 
hr! for both sub-runs Q.W.L.1.A and Q.W.L.1.B. 
and the heat transfer coeflicient value obtained 
was 342 B.t.u hr ft® F. 

However, in a few cases the values of heat 
transfer coeflicients varied considerably for the 
same set of conditions, because a minute variation 
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in the slope 6, of the line In At vs. 0, greatly 
affected the heat transfer coetlicient values, For 
example in run G.W.IL1 (Fig. 4) individual heat 
transfer coeflicient values were 560, 315 and 645 
B.t.u hr ft?/°F. These h,, 
averaged as they differed so greatly. It is however 
fortunate that the values of 6 had very slight 
differences given by 1558, 1519 and 1560 hr"! 
respectively. Instead of the 
transfer coetlicient which were differing 
appreciably, the averaging of the 6 values which 
varied slightly was carried out, and the heat 


values cannot be 


averaging heat 


values 


transfer coetlicient was calculated from the aver- 
aged value of the slope obtained in the triplicate 
runs for the same set of operating conditions. 
Thus in the above run G.W.II.1, the value of the 
155 B.t.u hr ft? F, 
from the averaged values of the slopes obtained 


coetlicient, was calculated 
in triplicate runs. 


The 


are presented in Tables 2 and 3. 


heat transfer coetlicients thus obtained 


It is seen from 














°F 





Temperature difference, 


¢~4. 











Fic. 3. Temperature difference vs. time (for ran Q.W.1.1) 
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these Tables that the heat transfer coetlicients 
between the solid particles and the fluidizing 
medium increased with the mass 


velocity, and in many cases they increased with 


increase in 


increase in particle diameter. 

The heat transfer coetlicient values obtained are 
also summarized in Table 4 and compared with 
the previous data on gas-solid heat transfer 
coellicients in fluidized beds. Values of Nusselt 
number (h, D, K) are plotted against Reynolds 
numbers (DG y) in Fig. 5 and the results are 
compared with Kerrenrine et al. [1]. It is found 
that the present values are much higher than the 
KeETTENRING data, because of the use of a liquid 
instead of a gas as the fluidizing medium. (Even 
though the “ over-estimated ~ heat transfer co- 
efficients [1], resulting from the assumption that 
the temperature of the gas was equal to the 


temperature of the wet fluidized solids, were used 
here.) The differences in both the results may 


also be due to the difference in fluidization 




















Temperature difference, 








(-*]. 























Fic. 4. Temperature difference vs. time (for run G.W.11,1) 
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Table 1. Observations of runs Q.W.AA.A and Q.W.LA.B.  Quartz-water systems. 
Quart: size; 0-00456 ft. 





QWLAA Q.W.1A.B 


Temp. reading Temp. reading 
Time on Beckmann ; on Beckmann 
(sec) thermometer > thermometer 
(c)* (cy 


OF 
1-4 
20 
23 
30 
$75 
+s 
50 
545 


1. Temperature—Time history 


Ss ss Ss OS 8S 
a ¢ = — a a 
a) oe 


ues = > 
— 


— 


Flow rate of water (Ib hr) 

Mass velocity (Ib hr ft*®) 15.030 
Pressure drop in main column (Ib) ft*) 23-5 

Pressure drop in glass column (Ib. ft®) “We 

Fluidized bed hieght in glass column (cm) 25-5 

«, porosity 0-68 
b, slope (hr-!) 2060 
h,, heat transfer coeflicient (Btu hr ft? F B42 





* The Beckmann thermometer was adjusted to read 0 C at 80 FP. 


Note > Zero of Fig. 3 corresponds to 20 sec. for run Q.W.1.1.A and 10 sec. for run Q.W.L.1.B. 


Table 2. Experimental values of heat transfer coefficents. Granite- water systems 





D, ‘ Modified Reynolds number h 
(ft) Porosity D, Gp (B.t.u hr ft® F) 


0-00360 0-680 10-7 
O-003560 O70 26-55 
O-00560 OT 204 
0-00300 708 16-8 
0-00300 727 210 
O-00250 728 1415 
O-00250 747 77 
0-00250 764 19-1 
O-O0212 747 We? 
o-oo2zl2 O77 4 
o-oo212 o7ol 17-8 
O-OO178 O-80 13-45 
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Experimental values of heat transfer coefficients. Quartz—-water systems 





€ Modified Reynolds number hy 
(ft) Porosity Dy Gp (B.t.u hr /ft®, F) 


0-00456 0-68 32-2 342 
0-00456 0-681 B4-7 409 
0-00456 0-702 38-5 550 
0-00456 0-710 $2-0 620 
0-00360 0-710 24-6 309 
0-00360 0-719 26-6 380 
0-00360 O-725 20-5 6 
0-00300 0-695 16-8 115 
0-00300 O-725 210 192 
0-00300 OTM 22-6 229 
0-00250 OT17 1415 129 
0-00250 727 17-7 385 
0-00250 0-747 19-1 474 
0-00250 0-766 23-1 S46 
0-oo212 O75 1? 145 
O-O0212 727 14-9 191 
o-oo2z12 0-748 16-08 141°. 
o-oo2z12 O70 17-8 503 
OOO1TTS O-850 12-5 183 
OMO1TTS 0-860 13-45 249 
OOO17T7TS O-874 14-7 304 
OOO1TTS O-R8S 16-3 309 





Table 4. Comparison of the present work with previous data on gas—solid heat transfer 


coefficients for fluidized beds 





Particles Range ofmass Max, range Effects of Effects of 
Type of system Technique used dia, range velocity of hy increasing increasing 
mesh size (Ib hr ft?) flow rate particle size 


Liquid fluidized beds Unsteady state 10-32 12,000- 20,000 113 620 Increases h,, Increases hy 
(Present work) 

Gas fluidized beds Unsteady state 16-115 100-600 0-06 0-9 No effect Increases hy 
(Wams ey et al. [4)}) 

Gas fluidized beds Steady state 14-18 0-750 3-0-1 Increases h, Increases hy 
(Kerrenrine et al. [1)}) 

Gas fluidized beds Steady state 28 48 5 920 33 69-6 Increases h,, Increases hp 
(Watton et al. [2}) 
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Fic. 5. Comparison of present data with 


KETTENRING ef al. (1). 

Mesh sizes Mesh sizes 

14-16 oO 10-14 

16-20 @ 14-16 

20-24 16-20 

25-28 10-24 

28-32 © 228 

OS ww-32 

mechanism. In Kerrenrine’s work the low heat 

transfer coefficients were due to the non-homo- 

geneous fluidization present in gaseous systems. 

The present values are higher than those of 

KeTTENRING because of the presence of homo- 

geneous fluidization, which eliminated the effects 

of the fluid by-passing, channelling and slugging. 

Also due to the absence of two phase fluidization 

the fluidizing liquid attains complete contact 

with the surfaces of the solids in the fluidized beds 

at any time. An equation similar to Kerrenrinc 
et al. [1] can be obtained from Fig. 5. 


», G\*? 
h D, @) “) 
fe 


D 
Pp"? — 0.00891 ( 
A 


which can be compared with the equation of 




















0,6 7 a 
“ ’ 
Fic. 6. Comparison of present data with 
Watton ef al. (2). 


Mesh sizes 
14-16 
16-20 
20-24 
25-28 
28-32 


Mesh sizes 
10-14 
14-16 
16-20 
10-24 
24-28 
28-32 


KetTTENRING et al. [1]: 
‘13 
Me Py _ vorss | Dy, “) (5) 
A im 


Equation (4) does not contain terms like voidage 
bed height, particle density, specific heat of the 
particle and other properties. 

The present data are also compared with the 
data of Waron et al. [2] in Fig. 6. There is some 
scatter in the results when (h, D,/K) is plotted 
against (D, G/y)'"(D,/D,)°*. However these 
values are less scattered than the data of WaLron 
et. al [2]. The equation for the present data is : 


h D,G 


Pp D 7 - O2 
— = 0-00571 “" (Dp/D,) (6) 


compared with equation of Warton et al. [2] : 
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heal tronsfer foctor 








Fic. 7. Relation between heat transfer factor and modified 
Reynolds number. 


Mesh sizes Mesh sizes 
V7 14-16 Oo Ww14 
V 16-20 @ 14-16 
6 20-24 16-20 
> 25-28 10-24 
O 28-32 © 24-28 


SO 28-32 


D, G 


hg D, 
kK Kt 


17 
(Dp/DJ°* (7) 


0-0028 


The low coefficients obtained by Watron may 
be due to (i) the averaging of the minimum and 
maximum limiting values of heat transfer co- 
eflicients taken and (ii) due to the by-passing of 
gas in the fluidized bed. Comparing Figs. 5 and 6 
it is observed that the results of equation (4) have 
less scatter than equation (6). 

The analogy between heat and mass transfer 
was also tried and it is found from Fig. 7 that all 
the present data are substantially below the 
theoretical curve given by Cuvu in Othmer’s book 
[7] and probably the relations given by Cau may 








heot tronsfer factor 








A» 








100 


Fic. 8. Comparison of present J, values with J, and J, 
values of KeETTENRING ef al. (1). 
Mesh sizes Mesh sizes 
JV 14-16 10-14 
V 16-20 14-16 
& 20-24 
> 25-28 
O 28-32 


not be satisfactory for liquid fluidized beds. While 
the present data is above the J, and Jp values of 
KETTENRING et al. [1] as presented by Cuvu in 
Orumer’s book [7], it is seen from Fig. 8 that the 
present results are nearer to the J, values of 
KETTENRING though their J, values were based on 
over-estimated heat transfer coetlicients. This was 
because of the high heat-transfer coeflicients 
attained in the present work. 


CONCLUSIONS 
The troubles encountered in the two phase 
fluidization are completely eliminated using a 
liquid as the fluidizing medium instead of a gas. 
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As expected, marked high heat transfer co- = Slope of the line In Af vs. 0 hr-? 


eflicients ranging between 113 and 620 = Heat capacity of liquid B.t.u/Ib/°F 
Heat capacity of solid B.t.u/ib/ oF 


B.t.u hr ft®, °F are obtained. The heat transfer Asvemmap Gensstnn of the yaatidle n 
cocllicients increased with increase in mass iininen af tien tain ft 
velocity and with increase in particle diameter. ;} = Mass velocity of the fluid Ib hr ft? 
Equation (4) seems to result in less scatter of Heat transfer coefficient between gas and 
solid particles B.t.u/hr ft? -F 
Heat transfer coeflicient between liquid and 
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KuLKARNI, Head of the Department of Chemical Tech- . Heat transfer factor dimensionless 
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the data than equation (6). 


Mass transfer factor dimesnionless 


ut a eats aie acai tate iaeeee abi Thermal conductivity of the fluid 
é “IS oO phice o COTRTL c sifieere app ‘ha- 
iuthors also sh phice PI B.t.u hr nr? F ft 
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Variable outlet liquid temperature 

Inlet liquid temperature 

Weight of solid particles 

Mass rate of liquid Ib /hr 
Viscosity of fluid Ib/hr/ft 
NOTATION = Void fraction in the bed dimensionless 
Time hr 


experimental set-up and work. 


Surface area of solid particles 
Temperature difference between outlet and ‘= ‘Time 
inlet fluid at the start of run P Difference in 


sec 
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Short-range production planning of continuous operations 


D. J. Witpe and A. Acrivos 


Department of Chemical Engineering, University of California, Berkeley, California 
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Abstract—A method is described for automatically scheduling production rates in certain 
systems of continuously operating manufacturing processes. The schedule which is suggested 
makes optimal use of availabie tankage in spire of random disturbances of supply and demand. 
Mathematically simple but exact, the method does not require knowledge of the statistical 
characteristics of the random fluctuations. 


Résumé—Les auteurs décrivent une méthode pour le prévision immeédiate des taux de produe- 
tion dans certains systémes de fabrication travaillant en continu. Cette méthode utilise au 
maximum le stockage disponible malgré les perturbations dues au hasard de la demande et de la 
fabrication. Methode mathématiquement simple mais exacte qui ne nécessite pas une connaissance 
des statistiques de fluctuations. 


Zusammenfassung—Es wird eine Methode beschrieben fiir die automatische Planung von 
Produktionsraten in gewissen Systemen von kontinuierlich arbeitenden Herstellungsprozessen. 
Die empfohlene Planung macht optimalen Gebrauch des verfiigbaren Lagerraumes bei zufalligen 


Stérungen von Angebot und Nachfrage. 


Die Methode ist mathematisch einfach, aber exakt und 


verlangt keine Kenntnis der statistischen Charakteristik zufalliger Schwankungen. 


INTRODUCTION 
Ix MANY chemical plants, between 10 and 20 per 
cent of the total capital expenditure is invested 
in tanks. The cost of the tanks themselves, 
together with the value of the property they 
occupy and the materials they contain, represents 
This investment 
produces no direct profit. It is in fact an economic 
liability, for it is continually subject to taxes, 


a sizeable outlay of capital. 


deterioration losses, and insurance charges. 

In spite of these drawbacks, tanks are a 
necessary part of any process, for they permit 
continuous operation of the system in the face of 
intermittent shipments, shutdowns, random dis- 
turbances and imperfect prediction of future 
supply and demand. As processing systems have 
grown more complex, production planning staffs 
have come into existence to handle the various 
scheduling problems which might be considered 
to fall into two categories—long range and short 
range. Long-range planners are concerned with 
setting average production rates over periods of 
time exceeding a month. In addition, a great 


deal of effort has recently been directed toward 
the solution of long-range problems in which the 
planner can focus his attention on optimizing 
such economic factors as profit or cost, ignoring 
those fluctuations in supply and demand which 
tend to cancel out in the long run. On the other 
hand, the short-range planner is left to figure out 
how to attain these long-range objectives in the 
face of day-to-day disturbances to the system. 
He is most concerned with the problem of avoiding 
situations which disrupt production, such as, for 
example, when some process has either empty 
feed tanks or full product tanks. From _ his 
standpoint, an optimal schedule is a flexible policy 
that can adjust to unforeseen random disturbances 
in such a way as to have as few major interrup- 
We shall be concerned here 
which 


tions as possible. 
with finding such strategies, we shall 


describe as optimal. 


Related to the operations problem described 
above is an interesting design question: how 
large should the tanks be? Too much tank 
capacity ties up capital, while too little leads to 
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many production upsets and greater need for 
Rather than attack this problem 
directly, we shall show how to find the minimum 


scheduling. 


tankage required when the short-range planner 
is given complete freedom to adjust production 
rates. This may be compared with the much 
larger capacity needed to absorb all fluctuations 
no adjustment rates is 


when of production 


permitted. By this comparison, the potential 

reductions in tank capacity through proper short- 

range planning can readily be appreciated, 
However, in spite of the potential savings in 


both capital and operating expense, little research 


has been directed toward the study of storage 
systems, possibly because an individual tank is 
such a simple piece of equipment that it has little 
intellectual appeal to an engineer. On the other 
hand, the analysis of systems of tanks is a challeng- 
ing theoretical problem, having practical interest 
to designers of new plants, to engineers engaged 
in process alerations, and to production planning 
staffs. 


PRELIMINARY SIMPLIFICATIONS 


A storage tank can have several functions. 
The major volume of the feed or product tanks 
must be reserved for holding material which is 
delivered or shipped intermittently but which 
must be processed continuously. Intermediate 
tanks are situated between processes to provide 


feed for 


upstream. 


downstream units during shutdowns 
The lowest part of a tank, where 
foreign matter collects, is not available for storage, 
as is the allowance for thermal expansion at the 
top. We shall be concerned with the volume re- 
maining after all of these allowances have been 
deducted. This buffer volume, as we shall call it, 
is available for absorbing the random disturbances 
in supply and demand rates external to the 
system. 

The production rates of actual processes are of 
course limited to a finite These limits 
complicate the scheduling problem considerably, 
and only certain simple limited rate situations 
have been solved [1]. For simplicity in this 
discussion it will be assumed that all production 


range. 


rates can be varied continuously over any range. 
The solution obtained under this strong hypo- 


1 


and A. 
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thesis will in many cases be a good first approxima- 
tion. 

The use of volumes, while helpful in visualizing 
the problems in storage systems, is inconvenient 
to work with because of density and stoichiometric 
In any 
system for which a complete material balance 


changes in the various processing steps. 


can be made, a given volume of substance at one 
point can be related directly to a weight of some 
other the To 
simplify the discussion, it will be assumed that 
all 


expressed consistently in weights of some reference 


material elsewhere in system. 


buffer volumes and tlow rates have been 
substance. With this done it is possible to picture 
each process as a pumping station transferring a 
homogeneous material, unchanged by processing, 


from one tank to another. 


‘ee - 
rur 


Consider a system of interconnected processes 


PROBLEM 


with many storage tanks, such as that shown in 
Fig. 4. 
prediction errors in the supply rates of raw 
As 


long as none of the tanks are completely full or 


There are random fluctuations and 


materials and in the demands of customers. 


empty, the system is considered to be operating 
properly. The problem is to find a policy for 
continuously adjusting the production rates of the 
that 
probability of satisfactory operation, in spite of 
the This 


paper derives such a strategy for the special case 


various processes will have the highest 


uncontrollable random disturbances. 
in which no physical restrictions are imposed on 
the system except that matter be conserved. In 
this simple but not unrealistic situation, the stat- 
istical character of the fluctuations need not be 
known. The optimal policy gives each instan- 
taneous production rate as a linear function of 
each random variable. 


Two 


The simplest non-trivial case arises when the 
system is a single process with a feed tank and 
a product tank, such as that shown in Fig. 1. 
Feed material from an outside supplier is stored 
in the feed tank. The process converts the feed 
into product which is stored for distribution to 
customers. Both the rate of supply and the rate 


SYSTEMS WITH TANKS 
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Supply rate uv 
Equivalents /unit time 








Production rate ¢ 
equivalents / unit time 


c- 


eae 
| Chemical 
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of demand are subject to unpredictable changes. 
It is desired to find out how to adjust the instan- 
taneous production rate in order to minimize the 
probability of filling or emptying either tank. 
Since it will be shown later on that any system 
can be reduced to this simple case, the two-tank 
system will be studied in detail. 

First 
instantaneous rates, average rates, and predicted 


we must make a distinction between 
rates. At any time @, let u denote the instan- 
taneous feed rate, in equivalent weight units per 
unit time, at which feed material is supplied to 
the feed tank. Over any time interval 0 — @, 
where @, is the starting point of the interval, the 
average feed rate 7 is given by 


° 


ud?’ (1) 


At time @,, at the start of a planning period, a 
prediction of i, the average feed rate at all future 
times 0 > @, is made. Let this predicted feed rate 
be denoted by @. 
while @ is a random variable. 


i is a non-negative constant, 
For abbreviation, 
we define the instantaneous feed rate prediction 
error Au and the average feed rate prediction error 
Ati by 

Au =u — i 

Ati 


= — 6 


In a similar manner let d represent the instan- 
taneous product rate at which product is shipped 
to the customer. d and d are defined as before, 
but we define the corresponding production errors 
slightly differently. 


l 


Customer 
demand rate 7 


Equivalents / unit time 





\ two-tank system. 


Ad 


Ad 


Thus the convention is adopted that a positive 
prediction error corresponds to a net gain in 
stored material over that predicted. 
Let f denote the instantaneous flow rate at which 
feed is converted to product in the process. f is 
defined as before. Since we have assumed that 
storage space has already been provided for all 
tank functions except buffering, we can specify 
without loss of generality that the predicted 


balances the predicted demand. 


i=-d=f 


The instantaneous and average flow rate deviations 


supply just 
That is, 


(2) 


are ; 
Ajmj—Js 


and 


f-f 


represents an increase in 


Af 
A positive deviation 
throughput. 

At time @ let W, be the amount of buffer 
material and E, the amount of empty space in 
the upstream (feed) tank, both in weight equiv- 
alents. The buffer capacity C, of the upstream 
tank is 


LE 


“u 


(3) 


Let W,° and E,° be the values of W, and E, at 
the beginning of the planning period (@ = @,). 
The changes in inventory AW, and in empty 
space AE, are given by 
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AW, = W, 
AE, = E, - 


Ww? (4) 
E,° (5) 
Similar quantities can be defined for the down- 


stream (product) tank and denoted by subscripts 


d instead of u. 
By taking material balances around each tank, 


we obtain 


f) da’ 


6 
AW, AE, | (u 
6% 


0,) (Ati — Af) 


AW, d) dé 


(0 — 04) (Af Ad) (7) 


We wish now toexpress Af as a function 4( Ai, Ad) 
of the random variables Af and Ad in such a 
way that the probability of satisfactory operation 
The 


operation are, of 


for every time @ > @, will be maximized. 


conditions for satisfactory 
course, that 
(8) 
and 
O< WA Cy (9) 


Any expression of Af as a function of Ai and 
Ad which maximizes this probability will be known 
as an optimal response, denoted by ¢*. 

It is clear now that any optimal response must 
be such that neither tank will go dry as long as 
there is any material in the other. It must also 
have the property that neither tank will fill up 
while there is any room in the other. We shall 
use these conditions to derive the optimal 
response. 

Let the average production rate at the time 
when either tank is empty be denoted by Af. 
Also, let Af, be the average production rate at the 
instant when either tank is full, Then, by com- 
bining equations (4) to (7) we obtain the following 
equivalent mathematical statement of the condi- 


tions which any optimal response must satisfy : 
¢* is optimal if and only if for all Afi and Ad. 


(A) Af, = Ati + W,°(@ — @,) if and only if 
Aj, = [Ad + W,°,(@ — @,)]. and 


and <A, 
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E,°, (0 — 0,4) if and only if 
E,, (0 — 0,)) 


(B) Af, = Aw 


Af; - [Ad 
When condition (A) is satisfied, 


(0 — 0) (At + Ad) + (W,2 + W,°) = 0 (10) 


Similarly, when condition (B) is satisfied, 


(0 — 0,) (Ati + Ad) —(E,° + E,°)=0 (11) 


The first term of equations (10) and (11) is the 
amount of material accumulated in the system 
at time 0. We shall call this quantity the aceumla- 


tion a. Thus, 


a=(0 04) (Ati Ad) (12) 


We shall require now that any optimal response 
be continuous in Ai and Ad, for any response 
involving discontinuous jumps in the throughput 
rate would probably not be too realistic and at 
any rate could always be approximated by a 
continuous function. This continuity assumption 
will lead to the rather surprising result that the 
response must be a linear function of the random 
variables, which naturally is very fortunate from 
the standpoint of simplicity. To prove this, we 
expand the optimal response function 4*( Ati, Ad), 
already assumed continuous, in a Taylor series 
about the origin, so that 


b* = ag + a, Ati + b, Ad + R, (Ad, Ati)... (18) 


where R, is the usual notation for the remainder, 
and a», a, and 6, are constants to be determined. 
Equation (2) implies that a, — 0, for when 
Ati = Ad = 0, 6* must also vanish. Furthermore, 
é* must satisfy condition (A). By equation (10), 

aa = — (aa 4 WP + We) 


14 
0 — 0, (14) 


when condition (A) is. satisfied. Combining 
equations (13), (14) and condition (A), we find 


that 


—(Ad+, 


ey 
+ b, Ad 4 R, (Ad. ad 5 


Since the left side of (15) is linear in both Ad and 
(WW? + Wy°)/(@ — O,)), the remainder R, must 


0 


ro 
Wy 


= —a(aa+ 


0 
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be identically zero. The coetlicients a, and b, are 


easily shown to be 


ay (16) 


Wes (17) 


°o 
W, 


and the 
therefore 


response satisfying condition (A) is 


~— 
Aj W,° Aw 


W,? 


W,® Ad 
=. (18) 

Ww, 

Since the response is linear in Ati and Ad, it also 

holds for the instantaneous variables Af, Au. and 

Ad. so that 

Wo? Au 

A P d 
J Ww, 


W® Ad (19) 

Ww? ? 
With this optional response then, the level 

changes are proportional to the accumulation, as 


can be seen by combining equations (6), (7), (12) 


) 


and (18), 

Wy? Ai 
WwW? 
Ad) 


,o 
AW, = (0 0,)| su W,? Ad 
(0 — 0) W,2 (Au 4 
Ww, Ww, 
aW,° 


AW. = ip i») 


‘ aw? 
AW w+ WwW, 
from which it follows that this optimal response 
distributes the accumulation in proportion to 
the amount of material the tanks at the 
beginning of the planning period. 

It of further to substitute 
equation (4) into equations (20) and (21) and solve 
for W,, and W,. One finds that 


from which, by dividing the first equation by 
the second, we obtain that 


is also interest 


"Ou ‘oO 
Ww, = We(i4 Set Me 
a 


°+ We 


W, = W? ( + Wi 


a 


1 


~ 
‘ 


at (22) 
Wa 

Therefore if this optimal response is used the 

quantities in the tanks are at all times proportional 

to the original amounts. 

It is apparent then from equation (22) that the 
linear response satisfying condition (A) satisfies 
condition (B) also only if by coincidence the 
original amounts are in the same proportion as 


the tank capacities, i.e., when 


wy Cc 
cc. 


= % 
Ww, 
In fact, we can show, by repeating our arguments 
used earlier, that the continuous response satisfy- 
ing condition (B), but not necessarily condition 
(A), is 
by? Ad 
ES 


Ej? Au 
Af BE» (23) 
Notice that the amounts W, and W, of condition 
(A) correspond to the empty spaces E, and E,, of 
condition (B). This duality makes it possible to 
use the results already obtained in equations 


(20), (21) and (22). 
AE 


“u 


(24) 
AE, 


and therefore, 
E 


E, 


E,° 


ae (26) 
E,? 


The difference between the two responses (19) 
and (23) proves that there is no analytic function 
satisfying both conditions (A) and (B). An 
analytic function, the kind familiar to most 
engineers, is one having all of its derivatives 
However, have not set any 
conditions on the derivatives of a _ response 
function, for we have merely required that the 
function itself be continuous. If we permit the 
first derivative to be discontinuous at only one 
point, we can find a family of continuous, non- 
analytic functions which satisfy both conditions. 
Such functions might be called piecewise linear, 
since the graph of the amount in any tank as a 


continuous. we 
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function of a, the accumulation (Fig. 2), is com- 
posed of two connected line segments of different 
slope. This graph must pass through three points : 
the initial condition and the two boundary 
conditions (A) and (B). Since the 
the break in the graph is arbitrary, a family of 
functions, rather than a unique function, will 


location of 


satisfy the conditions. 

Any optimal response causes the two tanks to 
behave as one large tank containing the total 
accumulation. takes ad- 
vantage of any cancellations in the fluctuations 
that may occur, as for example when both the 


Such an arrangement 


supply and the demand rate increase by the same 


amount. In an optimally controlled system, 
only excessive changes in the total accumulation 
a can cause unsatisfactory operation. However, 
since such circumstances are beyond the control 
of the planner, the use of an optimal response 
assures him that he is doing the best he can with 
the means at his disposal to prevent local emer- 
gencies. 

An example of an optimal piecewise linear 


response Is 


(W,? Au — W,° Ad) (WW? + W,%) 


when a 


,° @ 
E,°) 


when a 


(E,? Au — E,° Ad) (Eo 


This response is similar to the proportional control 
of servomechanism theory with a change in the 
sensitivity made every time the accumulation a 
changes sign. When the tanks are low, equation 
(27) distributes material; when they are high, it 
distributes empty space. 

An objection to the response (27) is that one 
might expect the shifting of constants to occur 
rather often, since for most probability distribu- 
tions of prediction errors, the original conditions 
have a high probability. Fortunately, the 
accumulation will actually tend to keep the same 
sign after the beginning of the planning period, 
This somewhat may 
apparent at first glance, and the sceptical reader 


surprising fact not be 
is referred to Feller’s second edition [2]. 

If both tanks are nearly empty to begin with, 
so that it is much more probable that they will 


amd A, 


\crivos 


run dry before they fill up, it is sensible to use 
the simpler response given by equation (19), even 
though it does not satisfy condition (B). A piece- 
wise linear variant of equation (19) is the following 


optimal response : 


CW? Au — W,° Ad 
Ww, — Wy 


when la| (WwW? — W,*) 


W,°) Au — (E,° —W,) Ad 
Ww? j ‘0 W,° 


“ul 


(Ww? 


ys 


(28) 
( kb? 
E,? 





ts ro 
i. when a W,,”) 


where again a is the accumulation defined by 
equation (12). 

Equation (28) is obtained by applying the re- 
W® Ad) (W,° + W,°) to the 
entire interval jal <(W,° + W,?) 
rather than, as in equation (27), only to the 
interval for which the accumulation a is negative. 
When a > (W,° — W,°) the proper coeflicients 
in equation (28) are as shown, since when 
a-~ W,? + W,? there are (E,° — W,°) units of 
empty space in the downstream half of the system 
and (£,° — W,°) units in the upstream half. 
Equation (28) is analytic in the symmetric interval 
ja| W° + W° and has the advantage that 
the probability of having ever to switch constants 
is quite low. However, if the tanks do fill up to 
the extent that condition (B) becomes important, 
the response will shift in time to avoid an overtlow 


sponse (W,° Au 
symmetric 


situation, 

On the other hand, when the tanks are almost 
full, empty space is critical and equation (23) 
is a reasonable response to use in spite of its 
failure to satisfy condition (A). An optimal 
piecewise linear response based on equation (23) 
is 
, EB? Au — Eo Ad 
EE, i E,° 

when ja] <(E,° + E,°) 
- (WLS — E,°) Ad 
Ww? 0 


(EB? 


- BE.) 


: ' (29) 
(Wy? — E,°) Au 
Wy — Ep 


when a< 





~ 


Equation (29) is dual to equation (28) and can 
be obtained from it by interchanging W,° with 
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~ Ef and W,° with E,°. The behaviour of 

the tank levels W, and W, as functions of the 
accumulation a are shown graphically in Figs. 2 
and 3 for the responses (23), (27), and (29) in the 
specific case where C, = 4, W,° » € 8, 
and W,° = 5. 














-6 a -2 0 
2, accumuiction, equivolents 


Fic. 2. Illustration of equation (27). 

















@, accumulation, equivolents 


Illustration of equations (23) and (29). 
equation (20) ------ equation (23). 


Fig. 3. 


LakGeE STORAGE SYSTEMS 

It will be shown now that any system of tanks 
and processes satisfying the simplifying assump- 
tions previously mentioned can be reduced to a 
two-tank system. This will make possible the 
extension of the two-tank theory already 
developed to large storage systems. 

In the discussion of large systems it will be 
convenient to use some concepts of graph theory 
[3]. The flow diagram for a system can be replaced 
by a graph in which the controlled flow rates 


are represented by directed lines (links) and the 
tanks by Streams with un- 
controllable or 


circles (nodes). 
rates 


into 


flow are 


the 


unpredictable 
indicated by dashed lines directed 
system. 

Figs. 4, 5 and 6 show how this representation 
is carried out. A conventional flow diagram for 
a hypothetical sulphuric acid company with two 
sulphur wells, two acid plants and three distribu- 
tion centres is shown in Fig. 4. Each plant has 
two units, the first to produce concentrated acid 
and the second to dilute some of the strong acid 
to make weak acid, One plant makes both grades; 
the other, only weak acid. There are two separate 
storage tanks at each distribution centre—one 
for each grade of acid. Demands for each grade 
at each storage centre and the supply rates of 
the sulphur from each well are random variables. 
Air and water are available as needed to make 
acid, 

In Fig. 5 the flow diagram is replaced by a 


Sulphur well 2 


Sulphur we’ 
+ 


Fieid2 

















Distribution 
centre | 








Fic. 4. A hypothetical acid company. 
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Simplified storage graph. 


graph. There are four points which one should 
especially notice: (1) The two sulphur storage 
tanks in plant B are combined into one equivalent 
tank because they have the same source and feed 
the same process. (2) Strong acid unit A is repre- 
sented by two lines because it has two sources 
of supply. The actual flow through the unit ts, 
of course, the sum of the flows through the two 
links. (3) Strong acid unit B and weak acid unit 
B are represented by one link because there is 














Fic. 6. Storage capacities and initial quantities. 
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no storage between them. They must both operate 
at the same equivalent rate. (4) The air and water 
streams are not indicated because they depend 
entirely on the acid rates in the units and hence 
are neither random nor controllable. 

Fig. 6 shows how the graph looks when it is 
untangled. Notice that there is only one loop, 
so that the analysis may be carried through 
without too much difficulty. 


Brancu Links anno Loop LINKs 


Consider now a link so located in a graph that 
if it were removed the graph would be separated 
into two sub-graphs entirely disconnected from 
each other. Such a link will be called a branch 
link. None of the links in a loop have this property 
because removal of such a link still leaves the re- 
maining graph in one piece. Such links will be 
called loop links. In Fig. 6, tank 6 is connected 
to tank 3 by a branch link and to tank 4 by a loop 
link. 

We shall now study large storage systems in a 
way identical to that used for the analysis of the 
two-tank problem. It is first immediately obvious 
that in large systems, as in two-tank systems, 
any response is optimum if and only if it guaran- 
tees that no tank will completely empty or fill 
before the others. It follows, therefore, that we 
may apply the two-tank results already derived 
to obtain the optimal response for any branch 
stream. This is so because all of the upstream 
tanks may be considered as a single large tank, 
since all of them must fill up or become empty 
simultaneously. This result also holds, of course, 
for the downstream tanks. In a similar manner 
the algebraic sum of the upstream prediction 
errors can be treated as a single random variable 


corresponding to Au in the two-tank case, while 
the analog to Ad is the algebraic sum of the 


downstream errors. Thus it then is a simple 
matter to find for each branch flow a piecewise 
linear response that will distribute material 
properly between the two halves of the system 
connected by the branch. It is of course necessary 
to specify such a response for every branch flow. 

As an example, we shall find the optimal 
response for the branch flow between tanks 3 and 
6. The equivalent two-tank system is shown in 
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Fig. 7, and the optimal piecewise linear response 
corresponding to equation (27) is 
(1390 
——— Aw, 
1610 
+ Ati 
+ Avy, 


220 
1610 
Au, - 
Aus) 
330 
t 2190 
Au,, + Atty, + 
Auy,) 


(Aug + Aty 


Auiys 
when a < 0, 

(30) 
1860 


—— Au 
2190 ° 


Au 


(Au, 


+ Aus 


when a > 0. 





a a Auy; 


Recall that a is the 
equal to the algebraic sum of the Au. 


accumulation, in this case 


4u,+ ds 

*JUp* Ju 

*JuUo* du: 
> 


Fic. 7. Equivalent two-tank system for link from ‘Tank 3 


to Tank 6. 


The problem of determining optimal responses 
for loop flows may be circumvented by converting 


the graph into one with no loops. called a tree. 
This is accomplished by arbitrarily removing 
enough loop links to break all loops in the graph. 
To “remove” a loop link all that needs to be 
done is to fix the flow rate in the link so that 
it does not vary. For example, the graph of Fig. 6 
can be converted into a tree by making constant 
the flow in the link from tank 8 to tank 13, 
Since this makes all of the links branch links, 
optimal responses for all of them can be deter- 
mined by the methods already discussed. Notice 
however that the optimal response in a loop link 
is not unique, for it depends on which link is 
removed to break the loop. 

Suppose, for example, we wish to find an optimal 
response for the loop link connecting tank 4 with 
tank 6. By fixing the flow in the link from 8 
to 13, we make the graph into a tree. With this 
link “ removed,” the upstream tanks are 2, 4, 5, 
8 and 14, making the required quantities 
W,° = 650, 2,,° = 900, WwW,’ = 960 and 

tq° = 1290. If instead it is decided to fix the flow 
in the loop link from 2 to 4, then only tank 4 is 


the 
hE? 


and quantities are 


1500 


upstream, 
WwW? = 110, 
E,f = 2080, 


necessary 


110, ~W,° and 


CONTROLLED 
SYSTEMS 


COMPARISON OF WITH 


UNCONTROLLED 


It is of interest to study next the behaviour of 
the hypothetical acid company with and without 
the optimal response. It may seem a little severe 
to compare a controlled system with one in 
which the flow rates are all held constant, but 
this latter practice is widely followed industrially 
because of its simplicity. The tanks are made 
oversize so that the unit rates can be fixed. 

Suppose that the prediction errors in the eight 
uncontrolled streams, indicated by dashed lines 
in Fig. 6, are each independently normally dis- 
tributed with a standard deviation of 10 weight 
equivalent units per day. The standard deviation 
of the average rate a is 
104/8 = 28-3 equivalents per day [1]. First we 
shall compare the probabilities of satisfactory 
operation for the controlled and uncontrolled 
systems on the thirtieth day after the beginning 


accumulation 


of the planning period. 

In the controlled system, the entire capacity 
of the system—-1610 equivalents of material and 
2190 of empty space—is the 
accumulation, which after thirty days has a 
850 equivalents. 


used to absorb 
standard deviation of 30 (28-3) 
Thus the controlled system is operating satis- 
factorily as long as the accumulation is neither 
1610 850 standard deviations below 
zero 2190,/850 above. By 
straightforward calculations, the 
probability of satisfactory operation can be found 
to be 0-966. 

On the other hand, in the uncontrolled system 
each tank must itself handle whatever fluctuations 
enter it. Tanks 3 through 8 will always have the 


less than 
por greater than 
probability 


same level, since they are not connected to any 
random disturbances, so only the remaining eight 
tanks limit the probability. After thirty days the 
standard deviation of each stream is (30) (10) 

300 units. Table 1 the individual 
probabilities of satisfactory operation for each 
of the remaining tanks. Since the fluctuations 
are statistically independent, the probability of 


shows 
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Table 


and A. 


ACRIVOs 


1. 





Tank 


Probability of satisfactory 


operation for the tank 0-500 


OM OsSsI O19 O-1S5 107 





satisfactory operation for the entire 
less than 1 per 


system Is 
given by the product of all these 
cent. 

Suppose now that we estimate how much 
tankage would be required in an uncontrolled 
system to give a probability of 96-6 per cent after 
thirty days. For simplicity, we assume that each 
outside tank is half full to begin with. Also all 
the outside tanks will be made the same size. 
The probability of satisfactory operation for each 
tank must be the eighth root of 0-966, which is 
0-996. This corresponds to 2-9 standard deviations 
on each side of zero, which means a total buffer 
capacity per tank of 2 (2-9) (30) 1700 equiv- 
For eight tanks, this is an aggregate of 
the 


almost to 


alents. 


13900 equivalents, compared to actual 


capacity of 2260. It is 
quadruple the capacity of the entire system to give 


necessary 


as good performance in an uncontrolled system as 
in an optimally controlled system. 


CONCLUSION 


We have shown how to obtain a set of simple 


scheduling rules which will maximize the prob- 
ability of successful operation for a storage system 
of any complexity, as long as there are no condi- 
tions to be met aside from tank capacity restric- 
The 
require a 
the 
The optimal scheduling 


tions and material balance considerations. 
of the 
the 

disturbances. 


construction rules does not 


knowledge of statistical character of 
external 
brings about one to one correspondance between 
the various tank levels and the total maount in 
the system, making feasible the useof feed back 
in the control system. The method proposed here 
makes use of the combined capacity of the many 
tanks in a system to absorb fluctuations, taking 
random 


In a_ particular 


advantage of any cancellations of 


deviations that might occur. 
example, it was found that optimal control gave 
a very high probability of satisfactory operation 
under circumstances in which conventional control 
would hardly work at all. A conventional system 
would need more than triple its actual capacity 
to give as good a performance as an optimally 


controlled system. 
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Abstract -Rates of heat transfer to non-Newtonian thuid systems passing through smooth round 
tubes under laminar flow conditions have been studied both experimentally and theoretically, 
for the case in which viscous heat generation may be neglected. The data of three independent 
investigators have been correlated with an average deviation of 16 per cent. The final equation 
contains the conventional Grashof number type of correction to account (empirically) for natural 
convection effects and the Sieder-Tate correction for radial variation of fluid consistency. In- 
terestingly, up to 37 per cent of the total heat was transferred by natural convection even in 


systems which were highly non-Newtonian, hence viscous, in nature. 


Résumé—Les auteurs font l'étude théorique et expérimentale des vitesses de transfert de 
chaleur de systéme fluides non Newtoniens s‘écoulant en régime laminaire dans des tubes lisses ; 
les données de trois auteurs différents montrent un écart moyen de 16°,. Dans l’équation finale, 
la correction apportée par le nombre de Grashof tient compte (empiriquement) des effets de 
convection naturelle, et la correction de Sieder-Tate de la variation radiale de viscosité, Cette 
équation est interessante jusqu’a 37°, de chaleur totale transportée par convection naturelle dans 
les systéme fortement non Newtoniens. 


Zusammenfassung—Die Wirmeiibertragung an nicht-Newtonsche Fliissigkeiten, die durch 
ein glattes rundes Rohr in laminarer Strémung fliessen, wurde experimentell und theoretisch 
untersucht, wobei Reibungswiirme vernachlissigt wurde. Die Werte von drei unabhingigen Fors- 
chern konnten mit einer mittleren Abweichung von 16°, dargestellt werden. Die Endgleichung 
enthilt die geowOhnliche Grashofzahl zur empirischen Beriicksichtigung der freien Konvektion 
und die Korrektiir nach Sieder und Tate fiir die radiale Verinderung des Fliessverhaltens. Intere- 
ssanterweise wurden bis 37°, der gesamten Wiirme durch natiirliche Konvektion iibertragen, 
sogar bei hochgradig nicht-Newtonschen und zugleich sehr zihen Systemen. 


Previous Work and THerory conditions (i.e. a well-developed, steady velocity 
Tuts paper deals with the rates of heat transfer profile) the presence or absence of visco-elastic 
effects is irrelevant as far as the velocity distribu- 
tion of the fluid, hence the heat transfer rates, 
are concerned [1]. Thus, while the analysis 
* purely- 
viscous,” the results are equally applicable to 
visco-elastic fluid systems until complications due 
to turbulence, natural convection or radial (and 
axial) variation of physical properties due to 


to non-Newtonian fluids flowing through smooth, 
round tubes under laminar conditions. The 


problem of interest is additionally defined by the 


following restrictions or assumptions : assumes the non-Newtonian fluid to be 


1. The fluid is not so viscous that heat genera- 
tion within the fluid need be considered. 
2. The physical properties of the fluid are 


invariant with respect to time, i.e. thixotropy is 
temperature changes are introduced, 


Table 1 reviews the recent studies in this area 
Under steady laminar (and isothermal) flow since about 1955; an earlier publication [2] has 


assumed to be absent. 


*Present Address: National Aeronautics and Space Administration, Cleveland, Ohio. 
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discussed prior contributions and incorporated 
any significant results thereof. All the theoretical 
approaches have assumed the physical properties 
of the fluid to be independent of temperature and 
have neglected contributions due to natural 
convection. The former assumption has already 
been shown to be a significant limitation [2] and 
the latter, as we shall see shortly, is responsible 
for the rather significant differences (10 and 25 
per cent) between our own earlier experimental 


Table 1. 


data [2] and those of Cuan [3] and Hirat [4], 
respectively. 

The over-all state of the art may be briefly 
summarized as follows : 

1. Paper number 3 in Table 1 serves to support 
theoretically the heat transfer coeflicients pre- 
dicted by paper number 4 in the only region of 
variables in which the latter were approximate 
rather than rigorous. Paper number 7 does not 
constitute any additional contribution as it is 


Review of previous work 





Major assumptions as to 


Authors i re - 
Heat flow 


. Hira [4} const, Bingham 


plastics (no assumption of a 
flow equation necessary in latter 


use). 


. Gricurt [5) const. 


3. Lycur and const, 


Birp [6] , 


const. 


. Merzner et al, [2] 


essary, provided No, 


Radial 


properties included empirically. 


. Cuarn [3] const. 


i}. SCHENK and const, 
Van Laar [8] law. 


WutreMan and const. Power law. 


Drake [9) 


. Wissier and const. 


Scnecrrer [10) 


. Brrp [11) const. Power law. 


plastics or 


Power law fluids. 


Power law fluids. 


General, no assumptions nec- 


variation of 


Power law or hyperbolic sine 


Bingham plastic fluid. 


Results 


Rheological properties 


pseudo- Duplicates earlier work of Picrorp 


[7] discussed in [2]. Has experi- 


mental data on slurries. 


Develops equations for point co- 
eflicients after shape of temper- 
ature profile is well established. 


Equations developed for both heat 
transfer rate and temperature pro- 
files, but only to Ne, = 82. 

Predictions rigorous if Ne, above 
20-100 or below 4, if n 0-10. 
other regions. 


20-100. 
physical Approximate in 


Checked experimentally. 
Work primarily experimental. 


Local Nusselt numbers computed 
for wide range of Graetz numbers. 
Heat generation included. 


Similar to [6] except for greater 
range of Graetz numbers.  In- 
accurate for highly non-Newtonian 


fluids. 


Computations available but not 
reduced for applicability to design 
problems. Heat generation in- 
cluded. 


Scope similar to [8] although heat 
generation not considered. 
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inaccurate for fluids having low flow-behaviour 
indexes and, at the higher Graetz numbers, 
merely duplicates the results obtainable directly 
from paper number 4 without the restriction of 
power-law behaviour. Of these, only paper 4 has 
considered heat transfer to dilatant fluids. 


2. Papers 1 and 5 provide some useful and 
badly needed experimental data (especially paper 
number 5) but do not add to the available 
theoretical studies. 


3. Papers 2 and 8 are of interest but of limited 
scope—the former because heat transfer in short 
tubes is not considered and the latter because 
the rheological equation chosen is, at best, an 
approximation to the behaviour of rather few 
real materials. 


4. Papers 6 and 9 cover the important area 
of constant heat flux calculations. Experimental 
confirmation is unavailable in this area, hence 
further work is obviously necessary. 


Returning to the problem of heat transfer with 
a constant wall temperature, it will be seen that 
corrections for natural convection must be in- 
cluded. This will be done empirically in view of 
the difficulty of any theoretical analysis directed 
toward solution of this problem when the tubes 
are horizontal [12], the configuration used to date 
in all experimental studies and probably in the 
majority of industrial installations. As a basis, 
the Merzner—Vavucun—Houcuton design equa- 
tions will be used, which may be summarized as 
follows when the effects of natural convection 
are excluded ; 


20 and n’ > 0-10: 


Os gt \i/s 
*) = 1.75 (“ C, 
Y A 


l. If Ng, > 


5 Bn’ rs 
4n 


where 


If capillary-tube viscometric data from which 
the flow behaviour index n' may be determined 
directly [13] are unavailable, 5 may be evaluated 
as the true shear rate at the wall, for the non- 
Newtonian fluid, divided by 8V/D. 


2. If Ng, < 20 and/or n’ < 0-10, then — 


O14 
re (22) Y = (3) 
kAM3 \ y; 

where ¢ represents the value of the product 
(h, D/k) (y¥¢/y)""* for Newtonian fluids at the 
same Graetz number. ¢ is given graphically in 
many references, including [2] and [12], and A'’* 
is given in reference [2]. 

Superposition of the effects of natural convec- 
tion will be carried out by assuming the correctness 
of the equation recommended for Newtonian 
fluids [14]* for Graetz numbers above 20: 


1-75 [N, 12-6 (Np, Ng, D/L)’ *}'? (4) 


sk 
By analogy, the effects of natural convection in 
non-Newtonian fluids perhaps may be incorpor- 
ated by adding a correction similar to the above to 
equation (1) or (3), giving, in the case of equation 


(**) 


1:75 (Ng, + 12-6 (Np. New D/L)?*}'”* (5) 


Ole 


The Prandtl and Grashof numbers contain vis- 
cosity terms evaluated at the wall temperature 
and shear rate since for most 
the viscosity will be lowest here, the resistance 


non-Newtonians 


offered to convective motion minimized and the 
over-all effect is greater than elsewhere in the 
tube. 
reduce convective motion, equation (5) must be 
considered to be limited to purely viscous fluids 
although equations (1) and (3) are not. Further, 
the viscous resistance to natural convection is in 
the transverse direction, rather than the longitu- 
dinal. To the extent that these may differ in non- 
Newtonian systems the correction chosen may be 


Since viscoelastic effects would tend to 


in error even in purely-viscous fluids. 


*A similar equation having the coeflicients 12-6 and 0-4 
replaced by 0-04 and 0-75, respectively, has been tentatively 
recommended by McApams [12]. This predicts coeflicients 
which differ appreciably from those of equation (4), 
however, and cannot be supported by the extensive data 
collected by EuBank and Proctor [14]. 


187 





\. B. Merener and D. F. Guock 





| 
Flow behaviour Diameter 


Symbol | luthor Fluid index of fluid | (in.) 


| GLUCK Carbopol O75 Poe 


@ 
A GLUCK ‘ arbopol on om 
° 


CHARM Nnunmonium alyinate Ow OOO 
| CHARM \pplesauce OOS 
CHARM Banana puree OM 


CHARM Banana puree OM 





E quaton (5) 


a 


ee = 








1000 


+126 (Me Move D/LI°* 
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EXPERIMENTAL CONFIRMATION AND 


Disc USSION 

The previously published data [2], on very 
viscous non-Newtonian systems, in which natural 
convection effects 
obeyed equation (1) over the range of Graetz 
numbers from 100 to 3,000 with a mean deviation 
of 13-5 per cent. Recent data of Cuarm [3] and 
of Giuck [15] on much less viscous fluids, fall 


must be negligibly small, 


appreciably above equation (1), as shown in 


Fig. 1, but correlate considerably better when the 
empirical for natural 
equation (5), is introduced (Fig. 2). The deviations 


correction convection, 
are indicated in Table 2. 


Table 2. Mean deviation, in per cent, from 





Equation (1) Equation (5) 
Data of CHarM 


Data of GuucKk 





The appreciable scatter, especially in Cuarm’s 
data, precludes any firm conclusions as the to 
actual correctness of the natural 
term chosen, although it is clear that some correc- 


convection 


tion is necessary and equation (5) is at least 
approximately valid. It is of interest, that if one 
were to delete arbitrarily the alginate data from 
Figs. 1 and 2 (which were Cuanrms’ earliest data, 
hence subject to the greatest random errors) the 
mean deviation of Cuarm’s data from Equation 
(5) wonld be increased to 13%, 

The only other available data are those of 
Hira [4]. These are on suspensions which are 
even less viscous than the fluids of Figs. 1 and 2, 
and the importance of the natural convection 
term would be increased. Qualitatively this is the 
case—most of the data fall 25 per cent above 
equation (1)—but they are not tabulated in 
sufficient detail to permit any more quantitative 
analysis. 

ConcLustions, DestGn RECOMMENDATIONS 

AND THE NEED FOR FurtHEeR Work 


Under conditions of constant wall temperature 
and laminar flow, heat transfer rates may be 


predicted within an accuracy of about + 10-15 
per cent, for viscous non-Newtonian fluids, by 
means of equations (1) and (3). The original 
development of equation (1) requires no assump- 
tion as to particular type of non-Newtonian 
This is not true in the 


region covered by equation (3), but in this area 


behaviour encountered. 


the possible range of the heat transfer coeflicients 
is so limited that an equal accuracy must neces- 
sarily be obtained. These equations presuppose 
the material is sufliciently fluid that viscous heat 
generation may be neglected yet not so fluid that 
natural convection effects become important. 

Natural convection may appreciably alter the 
predicted rates of heat transfer, even in non- 
Newtonian systems. Equation (5) correlates all 
the available non-Newtonian data and is recom- 
mended in this case. However, the data upon 
which it is based (Figs. 1 and 2) are not sufficiently 
extensive to indicate more than an approximate 
correctness for it. 

Under conditions of constant heat flux the equa- 
tions of Scuenk and Vaw Laar [8] and of Birp 
[11] are recommended as a first approximation. 
No data are available to confirm them and to 
the for 
variation of fluid consistency and for the effects 


develop corrections necessary radial 
of natural convection. 
As the for 


amounts to as much as 37 per cent and that for 


correction natural convection 
the radial variation of consistency 43 per cent 
(even for the limited data available, which were 
generally obtained under conditions in which 
these effects would be minimized) it is obvious 
that the need for further work is primarily in the 
direction of confirming or improving upon these 
terms. Conversely, it is suggested that the law 
of diminishing returns may be applicable to the 
mathematical analysis of oversimplified problems 
in this area since the effects of radial variation of 
physical properties and of natural convection may 
overshadow the importance of the non-Newtonian 
behaviour. 
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NOTATION 


(Any consistent set of units may be employed) 


Area of heat transfer surface Ratio of non-Newtonian to Newtonian fluid 
Specific heat of fluid wall shear rates, at a given flow rate 

Inside diameter of tube . Ratio of non-Newtonian to Newtonian heat 
Acceleration of gravity transfer coeflicients 

Heat transfer coefficient based upon an arith- Fluid viscosity. For non-Newtonian fluids, the 
metic mean temperature difference local ratio of shear stress to shear rate 

Fluid thermal conductivity Viscosity evaluated at wall conditions of shear 
Heated length of tube rate and temperature 

Flow behaviour index on’ = d (log ‘w! Fluid density 

Heat flow rate d (log 67/D) _ «= Shear stress at wall 

Wall temperature Grashof number. Ne, = 8 (At) D8p2g yp? 

Wall fluid (mean) temperature difference J Grashof number evaluated at wall conditions of 
Mean tluid velocity temperature and shear rate 

Mass flow rate N Graetz number. Ne, = wl, kL 

Volumetric coefficient of thermal expansion N Nusselt number. Ny = 4gD k 

Fluid consistency (see e.g. Refs. 1, 2, 13) } Prandtl number, Np, = Cpe & 

Fluid consistency evaluated at the wall temp- J Prandtl number evaluated at wall conditions of 
erature temperature and shear rate 
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Abstract—A simple tracer experiment carried out on a large underground water-pipe enabled the 
following to be determined: the volumetric rate of flow, the linear velocity in the pipe and the 
amount of silt which had settled-out in the pipe. The longitudinal mixing or dispersion coefficient 
was also determined. The experiment offers a novel illustration of the power of tracer methods 
to provide useful information about otherwise inaccessible regions. 


Résume—Une expérience avec traceur réalisée dans une conduite d'eau souterraine permet 
de déterminer le débit volumétrique, la vitesse linéaire, la quantité de dépdt et le coeflicient de 
mélange longitudinal. Cette expérience est une nouvelle illustration de l'avantage des méthodes 
avec traceurs pour obtenir des données sur des régions par ailleurs inaccessibles. 


Zusammenfassung—An einer grossen unterirdischen Wasserleitung wurde ein einfacher 
Tracer-Versuch zur Bestimmung folgender Gréssen unternommen: den Volumenstrom, die 
Lineargeschwindigkeit im Rohr und den Teil des Querschnitts, der sich zugesetzt hatte. Auch die 
Liingsdurchmischung oder der Dispersionskoeflicient wurden bestimmt. Der Versuch zeigt 
eine neue Anwendungsméglichkeit fiir Tracer-Methoden, um Information aus nicht zugiinglichen 


Bereichen zu erlangen, 


INTRODUCTION 


Tue meTHop of adding tracers to measure flow 
rates or “ dilution metering ” 
called, has for a number of years been used occa- 
sionally by civil engineers [1], [2]. More recently, 
there has been renewed interest in tracer experi- 
ments [3] as a means of investigating the mixing 
which occurs in a flow system. Mixing in one of the 


as it is sometimes 


most simple types of flow system, a long straight 
pipe, is important for both theoretical and prac- 
tical reasons, and previous experimental measure- 
ments have recently been correlated by LEVENSPIEL 
[4]. This correlation has shown that more ex- 
perimental work is required, particularly at high 
Reynolds numbers and in the transition region be- 
tween streamline and turbulent flow. 

This paper describes the application of tracer 
methods to a straight section of a large diameter 
water main for which the Reynolds number was 
greater than 10°. The function of this main is to 
return cooling water from an oil refinery to a 
nearby river estuary, and a fairly heavy deposit 


At the time of 
the tests, no provision existed for metering the 


of silt had been laid down in it. 


flow. The purpose of the tracer experiments was 
therefore threefold ; to find the rate of flow; to 
determine the depth of the silt deposit ; and to 
measure the dispersion coeflicient. The presence 
of the silt does detract from the significance of 
the dispersion coefficient calculated, but since 
data on straight pipes are so limited, it is con- 
sidered worth reporting. 

In the experiments to be described, the tracer 
was added at a point upstream and samples 
taken at both ends of the test section. Often in 
previous work [5] of this kind, emphasis has been 
placed on injecting the tracer rapidly as a compact 
slug. Although injection as a slug may simplify 
interpretation of the results, it is in fact necessary 
only when the concentration is measured at a 
single point down-stream. It is most undesirable 
that the success of this kind of experiment should 
have to rely on rapid injection, which is difficult 
to achieve with any certainty. 
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THEORY 


If a slug of tracer is introduced into a fluid 
flowing in a pipe, the distribution of tracer can 
be regarded as a wave which becomes increasingly 
dispersed as it passes down the pipe. It has been 
shown [6] that the longitudinal dispersion may 
be treated in the same way as diffusion and is 
determined by a dispersion coeflicient D, which 
in a uniform pipe is independent of position. In 
particular, the following equation in which the 
bulk movement of the fluid is super-imposed 
upon the dispersion, can be applied 


(1) 


The concentration ¢ is assumed to be uniform 
across a section of the pipe. 


Idealized tracer injection [7], [8] 
If the tracer is injected very rapidly such that 
in a 


at zero time the material is concentrated 


plane at 2 = 0, then the boundary conditions of 
equation (1) are 


0, ¢c =O for alla +0; 


x 


where Q 


Volume which tracer would occupy 
at unit concentration. 


A = Cross-sectional area of pipe. 

At any subsequent time ¢, the variation of ¢ 
with x follows a normal or Gaussian distribution, 
being symmetrical about the peak of the curve 
(a ut)* 

+ Dt 


@ exp 
24 v/ (2 Dt) ‘I 


The concentration of tracer as a function of 
time at a position 7 = /, for example a sample 
point, is very close to a normal distribution, 
provided D/ul< <1. Physically this condi- 
tion means that the concentration profile does 
not change appreciably in the time taken for the 


wave to pass the sample point. 


Non-ideal injection 


If the tracer is introduced over a finite period 


’.. Ler 


of time, then the concentration profile is broadened 
and may in addition be asymmetric about the 
peak. Sampling at a point downstream results 
in a curve of concentration as a function of time. 
Curves for two such stations, a distance L apart, 
provide the following information irrespective of 
whether or not the tracer is injected in the 
idealized manner. 

(i) Wolumetric rate of flow v. The 
tracer dM which crosses a section of the pipe in a 
time dt is the sum of the bulk flow and the flow 


mass of 


arising from longitudinal dispersion 


de 
dM ve. dt — DA = dt 


ou 


Substituting for de dv from equation (1) and 


integrating 


DA de dt 


u ot 


The last integral is ¢_ Cy Which is zero, and 
the 


zero [9]. Thus 


second integral also can be shown to be 


where O is the area under the ¢, ¢ curve. 

Summarizing, the volumetric rate of flow at a 
particular point is given by the amount of tracer 
added, divided by the area under the c, ¢ curve 
relating to that point. 

(ii) Linear velocity u. Uf the tracer distribution 
is symmetrical, then on the basis of the diffusional 
model, the position of the peak will travel with the 

Thus the linear 

= 
YT * given by 


t,) where ¢, and ¢, are the times at 


mean velocity of the fluid. 


velocity u, defined as u 


u=L/(t, 
which the peak passes the first and second stations 
respectively. 

Recently Aris [10] has examined the solution 
of equation (1) obtained by using the Laplace 
transformation. that whatever 
the form of the function representing the tracer 
distribution, the linear velocity can be found 
from the difference between the mean times at 


He has shown 
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which the wave enters and leaves the test section. 


In the dimensionless form derived by Arts 


(3) 


where /, and /, 


i.e. the linear velocity between two points is 
obtained from the mean times or the first moments 
of the areas under the two ce, ¢ curves. 

(iii) 
shown that likewise for any function, the disper- 


Dispersion coefficient D. Avis has further 


sion coefficient can be found from the difference 
between the variances of the concentration—time 
curves taken at the beginning and end of the test 
section. 

In dimensionless form 


(s,. D 


s)(")=2/r 29 


where S, and S, are the variances defined by 


i.e. the dispersion coeflicient between two points 
is obtained from the variances or second moments 
of the areas under the two c, ¢ curves. 


EXPERIMENTAL 

Choice of nickel tracer and method of analysis 

Cooling water used by the refinery is taken 
from the Thames estuary and its salt content 
varies with the state of the tide. It was anticipated 
that the maximum concentration of tracer at the 
second sample point would correspond to the 
amount added dispersed in about 10,000 ft® of 
water. This meant that if the amount of tracer 
introduced were to be kept within reasonable 


bounds, a substance was required which could be 
analysed to an accuracy of 0-1 p.p.m. in sea water. 
tried, the salt of 
fluorescein and nickel sulphate. found 
that fluorescein could be analysed to the required 


Two tracers were sodium 


It was 
accuracy using a Spekker absorptiometer, but 
preliminary trials indicated that some of the 
fluorescein added to the cooling water was pro- 
bably absorbed by small quantities of oil present. 
The nickel tracer was therefore adopted for 
quantitative work, 

A polarographic method was found to be the 
most satisfactory way of determining the nickel. 
Using an ammonium chloride-ammonia solution 
containing 0-1 per cent of gelatin as the supporting 
electrolyte, nickel concentrations of the order of 
5 p.p.m. could be measured with an accuracy of 

0-1 p.p.m. 
0-1 ml of 1 per cent gelatin added, and the volume 
made up to 10 ml with a 2N ammonium chloride 


5 ml of the sample were taken, 


2N ammonia solution. Part of this solution was 
then transferred to the cell of the polarograph, 
a crystal of sodium sulphite added to remove 
dissolved oxygen, and the current-voltage curve 
For nickel 

known diluted 
with sea water taken from the cooling water main, 


recorded, calibration purposes, a 


solution of concentration was 


and prepared for the polarograph by the same 


procedure. 


Distillation 
40 ft°/sec approx 


Ist sample ————- 
point 


Test section 
462 ft 





2nd sample ——+ 
point ’ 
(a) 


Fic. 1 (a). Arrangement of cooling water system. 
Tracer added at a. 


(b). Cross-section of main duct. 


193 





Arrangement of cooling water ducts 


Part of the cooling water system of a distilla- 
tion unit at the refinery of two 
open-topped cooling boxes in parallel (Fig. la). 
Water from each cooling box discharges over a 


consists 


weir into a vertical downcomer and thence into 
a closed horizontal pipe. This water is then 
joined by the remainder of the cooling water 
returned from the distillation unit, and passes 
into the main return duct which serves the whole 
From the weirs onwards, all the pipes run 
full. The cooling-water return main is made of 
concrete and is of octagonal cross-section (Fig. 1b), 
although the area available for flow is reduced by 
the deposit of silt. Downstream from the distilla- 


site. 


tion unit are two sample points situated 462 ft 
apart, spanning a straight uninterrupted section 
of the main, which runs about 18 ft below ground 
level. The side 


taken, are about 2 ft long and are accessible from 


arms from which samples were 
man-holes. 


Addition of tracers and sampling 

The tracers were introduced by rapidly pouring 
a solution into one of the downcomers of the 
cooling boxes immediately below the weir. This 
was the most convenient point, and had the 
advantage that the solution could be 
quickly, although there was obviously some back- 


added 


mixing due to the discharge from the weir. 

Fluorescein was used to give a visual indication 
of the time required for the nickel tracer to reach 
the sampling points. First a quantity of fluores- 
cein was dumped into the cooling box downcomer 
whilst water was being run continuously from 
the sample points. After a few minutes the water 
became green, the intensity rose to a maximum, 
and then diminished. The nickel solution was 
then introduced, and samples taken every 15 sec 
during the period of time indicated by the 
fluorescien test. The time taken to collect each 
sample was about 4 sec. 


RESULTS 
Concentration—-time curves 


In Fig. 2 the concentration of nickel in the 
samples is plotted against the time of sampling, 


. Ler 


measured from the instant at which the tracer 
was added. Both curves have broadly the same 
shape, but are noticeably asymmetric, the 
concentration rising quickly to a maximum and 
then decreasing more slowly. Although the 
tracer was added rapidly, it is unlikely that the 
idealized form of injection was achieved. The 





g/cm «10° 





Ni concentration, 























Time ofter injection, min 
(a) 
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Ni concentration, 














Time after injection, min 
(b) 


Concentration-time curves. 
(b) Second sample point. 


Fic. 2. 
(a) First sample point. 


tail on the curves probably arises from the flow 
pattern at the point where the tracer was intro- 
duced. Qualitatively it is consistent with the 
region immediately below the weir acting as a 
well-mixed vessel of comparatively small reten- 
tion time. Apart from the fact that the curves 
cannot be represented by an equation in analytical 
form, their asymmetric nature is no disadvantage. 
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Tracer measurements on a large diameter pipe 


Volumetric flow rate 

The flow rate at each of the two sample points 
was determined from the area under the con- 
centration—time curve using equation (2). These 
results are as follows : 


7,000 g NiSO, . 7H,O. 
0-96 g sec/1. 


Amount of nickel added 

Area under curve, Ist point 
0-88 g sec |. 
53-6 ft? /sec 
58-8 ft® /sec 


56-2 ft3 /sec 


Area under curve, 2nd point 
Calculated flow rate, Ist point 
Calculated flow rate, 2nd point 
Mean calculated flow rate 


Linear velocity 

The mean time of passage at each point was 
found from the experimental results themselves, 
/ dct / St, rather than from the curves, and 
(3) applied to the linear 


equation calculate 


velocity. 


5-4 min 
10-6 min 
462 ft 
1-48 ft ‘sec 


Mean time, Ist point 

Mean time, 2nd point 

Distance between sample points 
Linear velocity 


Dispersion coefficient 

Like the means, the variances of the 
curves were determined from the original 
numerical data according to conventional statis- 


two 


tical methods. The dispersion coefticient was then 
calculated from equation (4), with vo/V = uw/L 
and using the value of the linear velocity u 
calculated above. 

0-33 min? 
0-67 min® 
4-3 ft? /sec 


Variance, Ist point 
Variance, 2nd point 
Dispersion coefficient 


Amount of the silt deposit 


The amount of silt present can be calculated 
without necessarily assuming that it is uniformly 
distributed, as long as equation (3) may still be 
applied to find u. If at any point 


A = Free cross-sectional area available for flow, 


and 


A, = Cross-sectional area of empty pipe without 
deposit, then, 


Volume of silt in test section 
‘ e 
(A, — A). da A, . dx 
i, t 


but 4 .da = v.dtand v is constant by continuity 


Volume of silt 


Expressing the amount of silt as the volume 
fraction of the pipe occupied by the silt : 
vl 


Vu 


Volume fraction 
From the above results, taking the mean of 


38 ft? ; 
the geometry of the duct (Fig. 1b), 


V,/L = Me 


the volumetric flow rates, vu and from 


53-8 ft?. 


Volume fraction of silt 0-29. 


If it is assumed that the silt is uniformly 
distributed with a plane horizontal surface : 
Depth of silt z 2.9 ft. 


Discussion 

Finding a non-radioactive tracer which could 
be analysed to the required accuracy in oil- 
contaminated sea water posed a difficult problem. 
As far as the analytical techniques are concerned, 
nickel is certainly satisfactory, but there remains 
the question of whether any nickel was lost 
between the point of addition and the sampling 
points. The fluorescein proved unsuitable because 
such losses obviously did occur, probably by 
preferential absorption in traces of oil. Nickel 
might have been lost by precipitation, by exchange 
with iron from ferrous metal pipes and fittings, 
and by permeation into the silt deposit. If any 
loss did occur, the calculated volumetric flow rate 
would be higher than the true value, although 
the spread of the concentration—time curve, and 
hence the mean and variance, would be less 
affected. The calculated volumetric flow rate is 





in reasonable agreement with rough estimates 
of the cooling water used by the refinery, so 
any loss must have been small. At the second 
sampling point the calculated flow rate is about 
10 per cent greater than at the first, but this can 
hardly be regarded as significant in view of the 
accuracy of the measurements. 

The concentration-time curves, each with a 
tail, have the shape to be expected considering 
the backwash under the weir where the tracer 
of the 
curves, more frequent sampling, with each sample 


was added. To improve the accuracy 


collected over a shorter period, would have been 
but diflicult to 
owing to the inconvenience of handling sample 
The tlow 


rates at the first and second sample points differ 


desirable, achieve in practice 


bottles at the bottom of a manhole. 


by about 5 per cent from their mean value which, 
in view of the number of samples taken, is about 
the accuracy to be expected. 

the the 
dispersion coeflicient has been determined are not 


Obviously conditions under which 
ideal owing to the presence of silt and the un- 
certainty regarding its distribution, but never- 
theless it is interesting to compare the result with 
the recent correlation of Levensrres [4]. In this 
correlation, which strictly is for circular pipes, 
the dimensionless group D ud is plotted against 
Reynolds number to give a curve similar to the 
familiar friction factor vs Reynolds number plot. 
Assuming that the silt lies on the bottom of the 
duct with a plane surface, the equivalent diameter, 
defined by 4 free area perimeter, is 6-3 ft. 
Estimating the kinematic viscosity of the water 
as 0-75 es, the Reynolds number is 1-2 « 10°, 
On Levensrre.’s chart, reading from the full 
line which is the theoretical curve corresponding 


to the roughness factor for commercial pipe, 


Dud is 0-21 compared with the experimental 
value here of 0-46, Considering that the concrete 
duct will have a higher roughness factor than 
commercial pipe, and is not circular in cross- 
section, the comparison is reasonable. 

The results obtained in the foregoing experi- 
ments, while not of a high order of accuracy, 
provide a striking illustration of the usefulness 
and versatility of tracer measurements. The 
determination of the flow rate and the amount of 


. Ler 


deposited silt shows how tracers can provide 
information often 
without disturbing the normal operation of a 
plant. In this particular example, interruption 
of the cooling water system for inspection of the 


about inaccessible regions 


deposits would mean the whole refinery closing 
down. Only one section of the cooling water 
main was examined here, but by applying the 
same methods to other sections and repeating 
the measurements at suitable intervals, the ebb 
and flow of the silt deposits over a period of time 


could be followed. 
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NOTATION 
free area of cross-section of pipe 
area of cross-section of empty pipe 
concentration of tracer 
diameter or equivalent diameter of pipe 
dispersion coeflicient 


distance of sampling point downstream from point 
of tracer injection 


distance between sampling points 

area under concentration-time curve 

volume tracer would occupy at unit concentration 
variance of concentration-time curve 

time after injection of tracer 

mean of concentration time curve 

linear velocity 

volumetric rate of flow 

free volume of pipe between sampling points 
volume of empty pipe between sampling points 
distance downstream from point of tracer injection 


depth of silt 
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Abstract 
of solution of drops and bubbles under forced convection conditions. The effect of surface active 


Herein are reported experimental values and semi-theoretical formulae for the rates 


agents on the mass transfer rate is discussed and pertinent experimental data are presented. 
The following data obtained incidentally to the main study are also given: (1) velocity profiles 
in the entrance regions of pipes, and (2) diffusivities and other property value data for four 
water organic systems. 


Résumé 
les vitesses de dissolution de gouttes et de bulles dans des conditions de convection foreée. I 


L’auteur donne des valeurs expérimentales et des formules semi-théoriques pour 


discute l'effet dagents actifs de surface sur la vitesse de transfert de masse et présente des données 
expérimentales appropriées. I] donne également le profil de vitesse a lentrée des tuyaux ainsi 
que les diffusibilités et autres propriétés de quatre systémes eau-produit organique. 


Zusammenfassung—Fiir die Auflosungsgeschwindigkeit von Tropfen und Blasen bei erzwun- 
gener Konvektion werden Versuchswerte und halbtheoretische Formeln mitgeteilt. Die Wirkung 
von oberflichenaktiven Stoffen auf die Stoffiibertragung wird diskutiert undzu gehorige Versuch- 
sdaten angegeben. Ferner werden folgende Werte mitgeteilt, die neben den Hauptuntersuchungen 
erhalten wurden. 

1. Geschwindigkeitsprofile am Rohreinlauf. 

2. Diffusionskoeflizienten und andere Stoffwerte fiir vier Systeme von Wasser und organis- 

chen Stubstanzen. 


Many chemical engineering operations involve 
mass transfer between systems of spheres and 
surrounding fluid. A promising start towards 
a theory for mass transfer in such systems is 
the study of mass transfer between a_ single 
sphere and surrounding fluid. The results of the 
latter study may then be generalized to multi- 
sphere systems by additional work ; not even 
this will always be necessary since the rate 
expression for mass transfer from an isolated 
sphere occasionally well approximates the rate 
for a sphere in a multisphere system [1]. 

Mass transfer between a single sphere and 
surrounding fluid is usually divided into three 
parts : 


(1) Transfer within the sphere 
(2) Transfer through the interfacial phase and 
(8) Transfer outside the sphere. 


*Present address : 
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This paper attempts to characterize only a divis- 
ion of the last type of transfer, namely, forced 
The 
aim is to give reliable formulae for calculating 


convection mass transfer outside a sphere. 


the mass transfer coeflicient from such process 
variables as drop diameter and fluid velocity. 

The results of earlier studies directed at the 
same problem are summarized in Table 1, The 
present work significantly improves and extends 
the correlation of Nu, for spheres with surface 
flow and moderately improves the existing 
formulae for spheres with no surface flow. 

1. Tue ExXrPeriment 

The mass transfer coeflicient for the region 
outside a drop or bubble may be computed if one 
knows the average concentration of solute in the 
liquid outside the globule, the solute concentration 


American Cyanamid Company, Stamford, Connecticut, 
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at the interface and the rate of solute transfer. 
In the present work the solute concentration in 
the outside liquid was negligible. 

So that interfacial concentrations might be 
easily determined, the globule fluid and surround- 
ing liquid were chosen so as to form a partially 
soluble binary system with the fluid acting as 
solute in the liquid. If the usually reliable assump- 
tion of interfacial equilibrium is accepted, the 
the 
saturation concentrations for the binary system 


interfacial concentrations then equalled 
which are functions of temperature only. 
So that mass transfer rates might be easily 
measured, globules were formed and fed from 
fluid the 


diffusion of liquid into the globules was then 


saturated with outside liquid. As 
prevented, any change in globule size represented 
diffusion of solute away from the globule and 
coalescence of the associated saturating liquid 
with the outside liquid. Hence, the rate of solute 
transfer could be computed from the change in 
globule size or, equally well, from the amount of 
fluid fed to keep globule size unchanged. 

Much as in Ranz’s work [6] a drop or bubble 
was suspended on a glass capillary within a duct 
through which liquid flowed. As the drop or 
bubble dissolved in the passing fluid, its size 
was kept uniform by feeding fluid to it through 
the capillary from an assembly of ultramicro- 
burettes. The mass transfer was then obtained 
the 


measured intervals of time. 


from change in burette reading during 
In some runs fluid 
was not fed to the globule, but instead the mass 
transfer rate was found by measuring the rate of 
change of the globule diameter with time. Neither 
feeding of fluid nor change of globule diameter 
could have greatly perturbed the circulation 
within the drops or bubbles since the same transfer 
rates were observed in either case. 

The drops and bubbles were suspended slightly 
downstream from the entrance to a cylindrical 
duct. Upstream from the duct was a large tank 
in which the oncoming liquid was calmed. 
Adjoining the duct was a nozzle designed to 
give a uniform velocity over the initial cross- 
section of the duct. With this system the velocity 
field about the suspended drop was always 


laminar and nearly uniform in velocity. The 


velocity tields are described further in Appendix 1, 

Two methods were used to measure velocities 
in the ducts where drops and bubbles were 
suspended, With one method [21] the velocity 
field was made visible by colloidal tellurium 
trace particles and photographically recorded. 
With the other method [22 matter 
naturally present in the flowing liquid acted as 


colloidal 


trace particles, and these were viewed through 


an ultramicroscope. The preparations for and 
particulars of velocity measurement were essen- 
tially the same as those described by the originators 
of the techniques [21], [22]. Some experimenters, 
however, may find that the electrodes needed for 
the first mentioned method may be coated with 
tellurium more simply by electroplating, as done 
in the present work [23], than by the originally 
proposed vacuum evaporation technique, 

The mass transfer apparatus pictured in Fig. 1 
was used with bubbles 


for experiments gas 


dissolving in water. In experiments with drops 
dissolving in water it was unnecessary to de-gas 
the water; therefore the oxygen absorber and 
sand filter were deleted and the oxygen lines 
opened to the atmosphere. In experiments with 
water drops dissolving in isobutanol a similar 
but smaller arrangement than that pictured was 
constructed. 

Many of the physical properties of the chemical 
systems used were evaluated ; some of these data 
and the method of evaluation are reported in 


Appendix 2. 


2. Tunrory 
In this section the experimental system is 
idealized to a simple model. Equations for Nug, 
the dimensionless mass transfer coeflicient, are 
then presented for this model. The influence of 
surface active material on the mass transfer is 
also discussed. 


A. The model considered 


The essentials of the experimental system are 
embodied in that model where a solid or fluid 
globule is immersed in fluid and where the follow- 
ing conditions are fulfilled : 

(1) The interfacial tension is large enough to 
prevent oscillations and to keep the globule so 
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used with water as continuous phase and gas as 


Flow diagram of mass transfer apparatus 


dispersed phase 
A: Feed and collection tank. B: 
pump. C: Upper constant head tank. 
absorber. E : Sand filter. F : Rotameters. G : Calming 
tank. H: Upper flow nozzle and duct. I: Lower flow 
duct bubble. J: 
constant head tank. K : Dispersed phase feed system. 


Water circulation 
DD: Oxygen 


nozzle, and suspended Lower 


L. : Micam camera on microscope. M : Microscope lamp. 
N: Oxygen lines. O: Line to aspirator. 


nearly spherical that a single dimension, the 
diameter, characterizes its shape. 

(2) The velocity of the continuous phase is 
uniform in the region surrounding the globule 
except for perturbations caused by the presence 
of the globule. 

(3) The pressure and temperature are neither 
so high nor so low that the fluids cannot be 
treated as continua. 

(4) The fluids are so dilute in all components 
but one that: (a) binary diffusion equations are 
applicable, (b) heat generated by solution of 
transferring mass is small, (c) bulk flow in the 


radial direction at the globule surface is small, 
(d) the diffusivity, molar and mass densities, 
and viscosities do not vary significantly. 

(5) The pressure and temperature gradients 
are so small that they do not produce: (a) appre- 
ciable mass and energy fluxes, (b) appreciable 
variations in the physical properties. 

(6) Fields of forces other than gravity are so 
small that they do not affect the mass transfer. 

(7) The concentration of the diffusing com- 
ponents is uniform over the globule surface. 

(8) ‘Lhe process is steady state (this is infre- 
quent in commercial processes, but the actual 
unsteady state process may often be approximated 
by a series of properly chosen steady states). 

(9) There is no chemical reaction. 

(10) 


resistance to mass transfer and of negligible 


The interfacial phase is of negligible 


thickness (i.e. an interface). 

For the above model a dimensional analysis 
[23] of the applicable differential equations [24] 
and boundary conditions indicates that : 


Nu, = Nu, (Re,, Seg, 4/4. y*). 


Fortunately formulae derived for the above 
model can be applied to more general processes ; 
for example these formulae fitted the results 
from experiments where : 


(1) The ratio of the major to minor axis of 
the globule reached 1-4. 

(2) The concentration of solute reached 15 per 
cent by weight. 

(3) The ratio of the diffusivity near the globule 
surface to that at points far removed varied from 
0-5 to 1-5. 

(4) ‘Lhe bulk flow (the rate of solution of the 
globule divided by its surface area) approximated 
0-001 cm/sec, 

(5) The ratio of Gr,'/* ‘Re, reached 5. 


In the runs with large Gr, significant density 
differences existed. It has been suggested 
[25], [26] that a dependence of Nu, on Gr, should 
be noticed when the product Gr,'/* Re,'/* Se,"'’° 
exceeds 0-4; however, this constant is probably 
larger since no dependence was found in the 
present work even at the highest values en- 
countered, which were slightly above 1. 
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B. Surface effects 

The quantity y* enters the correlation when 
FRUMKIN and 
have suggested the most likely 


surface active matter is present. 
Levicn [27 
mechanism underlying the observed reduction 
in mass transfer rate by surface agents. Surface 
active matter accumulates on the downstream 
surface of a drop or bubble since on reaching the 
interface at upstream positions it is dragged by 
Then, 


because of the composition gradient created, an 


passing fluid to downstream locations. 


interfacial tension gradient also exists ; and the 
latter gradient opposes surface flow. 

The equation describing the hindering action of 
interfacial tension gradients to surface flow is a 
force balance at the interface. This balance is at 
any point on the interface : 


dy ox Tiy2 ' yrl (1) 


wherein y, X and y are respectively the interfacial 
tension, a vector distance parallel to the interface 


and the direction normal to the interface and 


Ty22 and 7,,, are viscous shears acting in the a 


direction at the interface and produced respec- 
tively by flow outside and within the drop or 
bubble. Equation (1) indicates that surface flow 
and, hence, circulation within a drop or bubble 
will be retarded appreciably and may even be 
0) when the 
tension gradient and external viscous drag are 


prevented (i.e. 7,,, interfacial 


of the same order. 


Derivations of Nu, 


The expressions for Nu, in several parts of the 
region where Nu, is a function of Re,, Sey, jy je, 
and y»* are presented in Table 2. The detailed 
derivations of these formulae presented 
elsewhere [23] ; only the general scheme by which 
formulae new in this paper are obtained will 


are 


be discussed here. 

Formulae 2 through 5 give Nu, for flows at Re, 
greater than 1 and with thin concentration 
boundary layers. They apply respectively to 
spheres with immobile surface, slowly moving 
surface, rapidly moving surface, and surface 
moving at potential flow velocity. 


Formula 2 is a correlation in terms of the 


pertinent dimensionless groups of the experi- 
mental data on rigid spheres taken during the 
work here reported. This formula differs from 
those advanced by others (see Table 1) primarily 
in that the large values of Se, in the present work 
showed the exponent of Se, to be nearer 0-35 than 
0-33. Formula 3 was obtained from formula 2 
on multiplying the forced convection term in the 
latter by a factor approximately correcting for 
slow interfacial flow. This factor was taken to be 
the ratio of mass transfer rates from plane plates 
with and without mobile surface ; the rate for a 
plate with mobile surface was obtained by the 
von Karman integral method as used by Eckert 
[28] for plates with immobile surface. Formula 4 
was obtained from formula 5 by replacing the 
potential flow interfacial velocity used in the latter 
by the actual interfacial velocity. 

Formulae 6 through 9 give Nu, for flows at 
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Fic. 2. The ratio £/ £, of the concentration boundary 
layer thickness with interfacial flow to that without 
interfacial flow. 
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' wet@TTteTS 


BY CONKIE AND SAVIC 


—-———— FROM INSERTING 
MEASURED VALUES 


OF NUp IN EQ'NS 3 AND 4 











V/ 
(ps,), REY? 


Fic. 3. The ratio kh, of the speed with which the interface between a sphere and surrounding fluid of real matter 


moves to that speed with which the interface would move if the sphere and fluid were frictionless, 


Re, less than 1. They apply respectively to 
spheres with rigid surface and thick concentration 
boundary layer (c.b.1.), rigid surface and thin 
c.b.1., slowly moving surface and thin c.b.1., and 
swiftly moving surface and thin c.b.1. 

Formulae 6 and 7 were derived by FriepLanDer 
[2] using the Stokes stream function [29] to 
describe the external flow. Formula 8 was derived 
using FrrepLaANpDER’s method but employing the 
HapaMarp stream function [30], which describes 
the external flow at Re, less than one about a 
sphere composed of any fluid and which reduces 
to the Stokes stream function if the sphere is solid. 
Formula 9 follows from equation (5) on replacing 
the potential flow interfacial velocity in the latter 
by the interfacial velocity predicted from the 
Hapamarp stream function. 

In order to use the mass transfer formulae 
derived for Re, greater than one, one needs 
expressions for: 


(1) € &, the ratio of the concentration 
boundary layer thickness with interfacial 
low to that without interfacial flow and 

(2) the ratio &, of the actual interfacial speed 
to the potential flow interfacial speed. 

The former ratio comes out of the derivation of 
formula 3 and is given in Fig. 2. Conkre and 
Savic [19] obtained the latter ratio for Re, 
greater than one from a balance between energy 
imparted to a sphere by the external flow and 
energy dissipation by viscous action within the 
sphere. Conkre and Savic’s values of k, are 
given in Fig. 3 together with the empirical line 
found by substituting experimental mass transfer 
data in formulae 3 and 4 and solving for k, ; 
the empirically found curve is recommended. 


D. Influence of y* on Nu, 
For reasons given earlier surface velocity is 
decreased where y* varies over the surface of a 
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direction of flow 


Y 


Fig. 4. 


hindered by interfacial tension gradients. 


Diagram of surface flow on spheres when 


fluid globule. Experimental observation suggests 
the action to be one of an immobile cap forming 
over the rear of the sphere while the surface flow 


on forward surface remains unimpeded (see Fig. 4). 
The extent of the cap increases as the variation 
of y* increases. Where the gradient in »* is 
sufficient to extend the immobile cap over the 
entire surface, the mass transfer relation becomes 
that for a rigid sphere. In this last case equation 
(1) may be integrated, since 7,,, equals zero and 
Tyee equals that for flow over a solid sphere, to 
give the limits above which rigid sphere behaviour 
obtains : 


Re, < 
Re, > 


Bo > 3 
Bo > 


1-36 Re,'/? 


3. EXPERIMENTAL Data 

The values of physical properties used in the 
dimensionless groups are arithmetic averages of 
the bulk and interfacial values. The interfacial 
values are for continuous phase fluid saturated 
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A graphical comparison between theoretically predicted and experimentally 


measured mass transfer rates for spheres free from interfacial tension gradients 
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with sphere phase fluid since equilibrium was 
assumed to exist at the interface and the sphere 
phase was initially saturated. 


A. Unhindered surface flow 


The pertinent experimental data are plotted in 
Fig. 5. The lines represent the values predicted 
from the formulae of Table 2 using the empirically 
determined curve of k,. The excellent fit of the 
data by the theoretical formulae is partly due to 
use of the empirical curve for k,; however, a 
very good fit is also obtained using the theoretical 
curve for k.. Further substance for the theories 
underlying the several formulae is given by the 
close approach at common boundaries of line 
segments representing different formulae. Formula 
4 is also in good agreement with data obtained by 
spray column techniques for water drops in 
isobutanol and the reverse [20], [31), for gas 
bubbles in water [32], and for the extraction from 
water of iodine or nitrophenol by carbon tetrach- 
loride [33). 


B. Zero surface flow 


The pertinent data are summarized in Table 3. 
In addition extensive data [5], [6] for drops 


evaporating in gases appear to have conclusively 
C Re,'* Se,", where 
C is between 0-55 and 0-60 and n could not be 


established that: Nu, — 2 


well established since Sc, was nearly one. Since 
the data of the first six systems in Table 3 are 
judged most reliable, these systems were used to 
compare formulae of the above type with n equal 
respectively to 0.33 and 0.35; a value of C 
between 0.55 and 0-60 was obtained only for the 
latter value of n thus establishing formulae 2. 
The remaining data in Table 3 are of uncertain 
reliability either because the diffusivity values 
are uncertain or because the gases involved 
(e.g. CO, and N,O) are so soluble that the experi- 
mental method was unsuited to the measurement. 

To obtain rigid sphere mass transfer data the 
globule surface flow had to be prevented ; this 
could be accomplished by adding oleic acid or 
powdered aluminum to the liquid drops. For 
low velocity flows past the drop and, therefore, 
for low shear stresses, surface active impurities 
already present in the liquids sufficed to prevent 


Table 3. The solution of rigid spheres of thirteen 

chemicals in water as described by empirical con- 

stants for the equations: Nu, = 2 + b Re,'* Se,'"* 
and Nu, = 2 + ¢ Re,'’* Se,°” 





Std. 
deviation 
of « 


Dispersed phase Runs 


ool 
ool 
ooo 
oot 
oor 
Nitrogen oor 
Air b 0-03 
O07 
O05 
oo 
OO 
oo 
oar 


Isobutanol 
Ethylacetate 
Cyclohexanol 
Isobutyl acetate 
Oxygen 


Carbon dioxide 
Nitrous oxide 
Argon 

Ethane 
Ethylene 
Hydrogen 





surface flow. With bubbles surface flow ceased 
soon after formation, presumably from adsorption 
of surface active matter present in the surrounding 
liquid, 


Partially hindered surface flow 


Examples of the mass transfer data obtained 
with partially retarded surface flow are presented 
in Figs. 5, 6, 7 and 8. The data of Fig. 6 were 
obtained by feeding the appropriate organic liquid 
plus oleic acid to a drop while maintaining a 
constant velocity flow of water past the drop. No 
quantitative measurement of Bo, the dimensionless 
interfacial tension difference, was made ; however. 
it is guessed that this plot is approximately one of 
mass transfer rate versus interfacial tension 
difference at fixed shear. 

The data of Fig. 7 were obtained by varying the 
flow velocity past a drop for which the hindrance 
to surface flow apparently stemmed entirely from 
surface active impurities in the chemicals used. 
An impurity depressing the surface tension as 
little as 0-1 dyne/em would completely suppress 
surface flow by the criterion given at the close of 
Section II, namely, Bo > 1-36 Re,'/* ; hence an 
even smaller impurity would suffice for the 
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Fic. 6. Reductions in the solution rates of drops by additions of oleic acid. 
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Fic. 7. Solution rates of ethyl acetate drops in water when hindered by interfacial tension gradients. 


207 





R. M. Grairrira 








6-7/e-—9] 10-1) 
(CM/SEC)| (CM/SEC)] (CM/SEC) 
9° o- 6 
V > v 





























¢ 
Cat BUBBLES, EQ'N 2, NO 


— 





o 


SURFACE FLOW 


- 





| 





6 


t/Voq (SEC*/CM) 


Fig. 8. 
behaviour shown in Fig. 7. Since the concentra- 
tion of surface agent in the chemicals was pre- 
sumably constant during the data collection, the 
interfacial tension difference is guessed to approxi- 
mate a fixed value for the data of this figure. 
Hence Fig. 7 may approach a plot of mass transfer 
rate versus viscous shear at fixed interfacial 
tension difference. 

In Fig. 8 the mass transfer rate of an oxygen 
bubble in distilled water (de-gassed by passage 
through de-greased steel wool) is given as a 
The 
rate is presumed to fall off as surface active 


It 


may be noted that the transition from essentially 


function of time after bubble formation. 


matter is adsorbed from surrounding water. 


fluid sphere to rigid sphere behaviour required 
about one minute; since this is the order of 
magnitude of residence times for bubbles in 
various equipment, similar phenomena may occur 
there. 


4. GuIpDEsS For tHe Use or TaBLe 2 


The equations given in Table 2 may be used 
for systems where drops or bubbles have either 


Reduction in the solution rates of oxvgen bubbles by gradual adsorption of surface active matter. 


zero surface flow or unhindered surface flow. 
If one suspects that surface agents are present 
but in amounts insuflicient to prevent surface 
flow, he may still use the equations of Table 2 to 
obtain upper and lower bounds on the mass 
transfer rate. One may note that rapid globule 
rise or fall, low interfacial tension, short residence 
times, and pure materials all reduce the likelihood 


of significant surface active effects. 
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APPENDIX 1 


Mass transfer studies and velocity measurements were 
carried out in the entrance regions of pipes with circular 
cross-section, Fig. 9 shows the velocity pattern achieved 
when the container emptied by the pipe was essentially 
infinite compared to the pipe diameter. Fig. 10 shows a 
velocity pattern achieved when the pipe was preceded 
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. Dimensionless plot of velocities in the entrance region of a duct like that shown in the legend. 
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by a relatively small container (i.e. a section of larger 
diameter pipe). Fig. 10 is very nearly the same as Fig. 9 
if Z 4 is substituted for Z in the latter. The quantity 4 
is presumably the fraction of the large diameter pipe length 
Z(=—lSem, Z d= 2-88) which added to 
Z when the ratio of d’/d, therefore, of large pipe diameter 
over small, is 1-40. The lines obtained by Nikurapse [34] 
from his experimental work are plotted in Fig. 0 and agree 


must be 


well with the present data. 

All sections of pipe had bell-shaped entrance nozzles. 
Several layers of I4-mesh screen were placed over the 
entrance to the large diameter pipe or over the small 


diameter pipe if the former pipe was not used. The data 


Table 4. 


in Fig. 9 were obtained with isobutanol in the flow system 
and using the ultramicroscope technique for velocity 
The data in Fig. 10 were obtained with 
system and using both the ultra- 


measurement. 
water in the flow 
microscope and tellurium techniques of velocity measure- 


APPENDIX 2 


The data reported in this section were obtained by 
The following tables identify the 
general technique used, present the data or equations 
determined by them, and give the estimated accuracy of 


well-known methods. 


these data and equations. 


Identification of the experimental methods 





Property 
Density 
Diffusivity 


Derivative of diffusivity with respect 


to concentration 
Interfacial tension 
Mutual solubility 


Refractive index 


Viscosity 


Measured with or by: 


Ravleigh and or Gouy interferometry 


Ravleigh interferometry 


Estimated accuracy (%) 


Ostwald Sprengel type pyenometer O1 


O54 


See Creern [35] 


Drop weight apparatus 5 
\nalvsis of co-existing phases in equilibrium 1 


Immersion refractometer 


Ostwald type capillary viscosimeter 





Table 5. 


p = A — B(25 


The density of several chemical systems as described by empirical constants for the equation : 


. T) C p4/p 





Solute Solvent 
Water 

Water 

Water 

Water 
Isobutanol 
Ethyl acetate 
Cyclohexanol 
Isobutyl acetote 


Isobutanol 
Ethyl acetate 
Cyclohexanol 
Isobutyl acetate 
Water 

Water 

Water 

Water 


O-9970 
O-0970 
O-9970 
9970 
(7980 
O-S8952 
00-9446 
O-S8665 





*In Table 5 p is the density of solution in g ml; T is the temperature in °C, 


and p,4 is the concentration of solute in 


2/100 ml. The constants have been checked only from 20 °C to 35 C and from p 4 equal 0 to its value at saturation, 
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Table 6. Diffusivity and the first derivative of diffusivity with respect to concentration for several chemical 
systems* 





A (Solute) B (Solvent) Dap 1,6 
Water 6-800 
Water 218 
Water 211 
Water “SRS 
Water 245 
Water 659 
Water TOl 
Isobutanol 2-068 
2-310 


Isobutanol 
Isobutanol 
Isobutanol 
Ethyl acetate 
Cyclohexanol 
Isobutyl acetate 
Isobutyl acetate 
Water 


Water Tsobutanol 


Dap #4 * PA 


Method 


O85 
5-92 
5-92 
5-92 
2-863 
0-298 
0-302 
12-26 


10-80 


G 
R 
G 


G 
G 


www wee heh Oe 
St Se Se Se Oe Se Se Se Ge 





*In Table 6, D is the diffusivity in cm?® see ; 


p4 is the volume weighted average concentration of the diffusing com- 


ponent 4 in the scluticn in g 100 ml; T is the temperature in C ; and G and R stand for measurements with the Gouy 


and Rayleigh interferometers respectively. 


Table 7. 


Interfacial tension at 25°C for four 


chemical systems 





System Interfacial tension (dyn em) 
2-132 
i410 
O85 


13-72 


Isobutanol water 
Ethyl acetate water 
Cyclohexanol water 

Isobutyl acetate. water 


Table 9. 
in water as described by empirical constants for the 


pP4 Pp A BR 


The refractive indices of three chemicals 


equation™ ; 





Solute Solvent A x 103 B 1 
0-962 
1-32 


0-709 


Water 
Water 
Water 


Tsobutanol 1-28202 
175918 


0-94483 


Cyclohexanol 
Ethyl acetate 





Table 8. 
as described by empirical constants for the equation® ; 


pa/p = A — B(T — 25) 


Mutual solubility of four chemical systems 





Solute Solvent 

8-11 

745 

3-85 36 
0-666 O79 
16-04 76 
3-33 240 
11-44 21 


1-22 


Water 
Water 
Water 
Water 
Tsobutanol 


Isobutanol 
Ethyl acetate 
Cyclohexanol 
Isobutyl acetate 
Water 

Water 

Water 

Water 


Ethyl acetate 
Cyclohexanol 
Isobutyl acetate 





*In Table 8 T is the temperature in C; py is the 
concentration of solute in the saturated solution in 
g/100 ml; and p is the density of the solution in g/ml. 
The constants have been checked only from 20 C to 


35 CC. 


*In Table 9 p 4 is the concentration of solute in solution 
in g/100 ml; p is the density of solution in g/ml; and 
R is the refractive index at 25-00 °C. The constants were 
checked only for p4 from about 2 units below saturation 
to saturation. 


Table 10. The viscosity at four temperatures of 
several chemical systems as described by empirical 


constants for the equation*: = A+ Bp, p 





A B 
1-008 
0-896 4-34 
0-801 3-64 


0-722 2-95 


10? | T 

19-94 
24-92 
29-96 
35-04 


Tsobutanol Water 


5°25 





of the solution in 
centipoise, p is the density of solution in g/ml, p4 is the 
concentration of solute in g/100 ml, and T is the tempera- 
ture in “C at which the constants are valid. 


*In Table 10 y is the viscosity 
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Table 10—continued 





Solute Solvent 
Ethylacetate) Water 1-008 288 19-04 
and isobutyl 0-896 2-: 24-92 
acetate 0-801 S! 20-06 
O-722 ". S544 
Cyclohexanol Water 1-008 54m 10-4 
0-896 : 24r 
0-801 33-3 20-06 
0-722 276 35-04 
Watert Tsobutanol 3-205 3-3: 19-94 
2-306 276 24 
2-560 “4 2068 
2-209 ¢ B54 
Water Ethylacetate 0-449 19-94 
| 0-424 76 24-92 
0-399 7: zoo 
O-378 S546 
Watert Cyclohexanol = 42-58 19-4 
31-47 24-0r 
23-72 29-06 
18-21 ' 35-04 
Water Isobuty! 0-691 7 19-4 
acetate 0-649 : 24 
0-609 20-06 
0-570 35-04 





+The constants for this system were checked only for 
p4 from 14 to 18 g/100 ml, 

tThe constants for this system were checked only for 
p4 from 8 to 12 ¢/ 100 ml. 


NoTaTIioN 


Bond number defined as 
y* (0 = 2) — y* (@ = 0) 
heat capacity per unit mass at constant pressure 
diameter of sphere (except in Appendix 1 and 
Figs. 9 and 10 where d is the diameter of the 
flow duct) 
binary diffusivity defined, for example, for com- 
ponent A in phase 2 as D 4, 
Grashof number defined as d® g p(p, pi) 
acceleration of gravity 
coefficient of thermal conductivity 
mass transfer coefficient defined, for example, 
in phase 2 as: Rag = N4/(P 420i — Paw) 

= dimensionless quantity defined as 
2v,/(3v,.. sin ) 
mass flux of species A based on a unit of sphere 
surface 
Nusselt number based on an area weighted 
average of the transfer coeflicient over the 
particle surface: Nu = kd/D 


Subscripts 


A 
i 


“ 


Peclet number for mass transfer defined as the 
product Re Se 

Prandt!] number defined as pep /k 

in Fig. 6 the total volume of liquid fed to the 
drop including that for formation 

Reynolds number defined as do, v (in Figs. 9 
and 10, » is evaluated for the fluid in the duct 
and ov, is replaced by V) 

in Figs. 9 and 10, distance measured in 
evlindrical co-ordinates from the axis of the 
duct 

dimensionless quantity defined as 2rd 

Schmidt number defined as » D 

in Fig. 8 the time after bubble formation 

local mass average speed 

an average speed defined as the total volumetric 
flow rate through a duct divided by the duct 
cross-sectional area 

in Figs. 9 and 10 and Appendix 1, distance 
measured from the end of the flow nozzle along 
the axis of the lower flow duct. 

interfacial tension 
quantity defined as y jg ge, 

the angle in spherical co-ordinates measuring as 
shown in Fig. 4 

the coefficient of shear viscosity 

kinematic viscosity (jy /p) 

quantity defined as the concentration boundary 
layer thickness divided by the hydrodynamic 
boundary liver thickness 

14150 

mass density of pure component (g of A/ml 
of pure 1) 

mass density of component .1 in solution 

(g of A mil solution) 


a chemical species 

quantity evaluated at the interface 

(1) quantity evaluated for spheres with immobile 
interface (2) or when applied to @ it designates 
the lowest value of @ at which surface is immobile 
quantity evaluated in (or for) the phase compos- 
ing the sphere 

Quantity evaluated in (or for) the continuous 
phase, where the latter is that phase surrounding 
the sphere 

quantity evaluated in phase 2 at a point 
sufliciently far from the sphere as to be un- 
affected by the sphere 


Superseripts 


quantity is dimensionless 


Marks above symbols 


quantity is a mass average value if a bulk 
property and an area average value if a surface 
property 
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Transfer from a sphere into a fluid in laminar flow 


M. Linron* and K. L. Suruertanp*® 


Division of Physical Chemistry, C.S.1.R.0. Chemical Research Laboratories, Melbourne 
(Reecived 23 Neowember 195%) 


Abstract——Boundary layer theory and experimental results for the transfer of heat and mass 
from spheres in forced flow are compared over the range of conditions : 


1 < Re Is 1 OS rr we 2 Nu 5M). 


The overall and local solution rates of 3.8 in. diameter spheres of benzoic acid were measured 
in uniform flow in a water tunnel! at 400 Re 7580. The overall transfer was correlated by 


Nu | 582 Rel? Prl/4 


The relative distribution of local transfer rates over the front half of the sphere was in fair agree- 
ment with theory and other workers but the absolute local values at the front stegnation pot 
were about 40 per cent lower than theory and differed considerably among other workers. The 
discrepancies may in part be due to a gradual transition from viscous flow to potential flow 
outside the boundary layer as Re increases from 1 to 105 The relative transfer from the rear 
half of the sphere was observed to increase with Re, The angle of separation 0, found experi- 


mentally varies approximately according to 
0, = 83 + 191 Re-¥/3 


The average transfer rate over the front of the sphere is theoretically proportional to 
Re’? Pr¥3_ Over the whole sphere 0-5 < Nu Re’? Pr-¥3 < 10 for 10? < Re < 10°, The 
approximation is only rough because of the factors dependent on Re mentioned above, and the 
following secondary factors: a lower limit to transfer rate set by diffusion or conduction, 
density convection unless Gr Re®, turbulence and or non-uniformity of the free stream, 
interference from sphere supports, influence of duct walls, uncertainty in transfer properties 
of fluids, variation of fluid properties with concentration or temperature, and dust accumulation 
on the sphere. 


Résumeé-——La couche limite théorique ct les résultats expérimentaux pour le transfert de chaleur 
et le transfert de masse a partie de sphéres dans un écoulement forcé, sont comparés dans 
Vintervalle : 

1< Re < 18 « 1° +05 < Pr< 108 + 2 < Nu < 5, 


Les vitesses de dissolution locale et globale de sphéres d'acide benzoique de 3 Sin. de diametre 
ont été mesurées dans un tunnel d'eau A 490 < Re < 7580. Le transfert global est donné par 


Nu = 0.582 Re! ? pri 4 


La distribution des vitesses de transfert sur la moitié avant de la sphére s'accorde bien avec 
la théorie et les autres expérimentateurs, mais les valeurs locales absolues au point de stagnation 
avant, sont environ 40°, inférieures 4 la théorie et varient avec les expérimentateurs. Les 
différences sont dues en partie 4 la transition graduelle de Pécoulement visqueux en écoulement 
potentiel en dehors de la couche limite quand Re augmente de 1 a 10°. 

Le transfert relatif sur la moitié arri¢re de la sphére augmente avec Re. Le point de séparation 
@, trouvé expérimentalement varie comme 0, = 83 + 191 Re~?®, 

La vitesse de transfert moyenne sur le front de la sphére est théoriquement proportionnelle 
& Re? Pr? Sur la sphére entiére 0.5 < Nu Rew? Pr? < 1,0 pour 10? < Re < 10°. 
L'approximation est un peu grossiére & cause des facteurs dépendants de Re déja mentionnés 
et des facteurs secondaires suivants: limite inférieure de vitesse de transfert par diffusion ou 
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Re?) turbulence et non uniformité 
influence des parois de la conduite ; 


conduction ; importance de la convection (sauf pour Gr 
du courant libre ; interférence des supports de sphéres ; 
incertitude des propriétés de transfert des fluides, variation des propriétes du fluide avec la 
concentration ou la température, et accumulation des poussiéres sur la sphére. 


Zusammenfassung—Fiir den Wiarme— und Stoffiibergang an Kugeln in erzwungener Stromung 
wurden die Ergebnisse der Grenzschichttheorie und der Versuche miteinander verglichen und 
zwar im Bereich von 1 < Re 18 « 10°; 05 < Pr < 108%; 2 < Nu < 540. 

Die gesamte und die Ortliche Auflésung von Kugeln aus Benzoesiure( Dmr. 9,5 mm) in 
einer gleichférmigen WasserstrOmung bei 490 < Pe < 7580 wurden gemessen. Die Gesamtwerte 
der Chertragung konnten durch folgende Bezichung wiedergegeben werden : 


Nu = 0.582 Re! ? Pr’ 3 


Die relative Verteilung der Ortlichen Ubergangskoeflizienten auf der Stirmseite der Kugel war 
in guter Cbereinstimmung mit der Theorie und anderen Versuchsergebnissen, aber die absoluten 
Ortlichen Werte am vorderen Staupunkt waren etwa 40°, unter denen der Theorie und differierten 
betriichtlich mit fremden Versuchsergebnissen. Die Unstimmigkeiten mégen zum Teil durch 
einen schrittweisen Cbhergang vom zihen Fliessen zur PotentialistrOmung ausserhalb der Grenz- 
schicht beim Anwachsen von Re von 1 bis 10° verursacht sein. Die relative Cbertragung an der 
Riickseite der Kugel stieg mit Re an. Der Winkel des Ablésungspunktes 6, dinderte sich nach 
den Versuchsergebnissen ungefilir nach der Gleichung 


“0 83 + 191 Re-!3 


Die mittlere Ubertragung an der Vorderseite der Kugel ist nach der Theorie proportional zu 
Re! 2 Pr! 3. Cher die ganze Kugel ist 0,5 < Nu Rew? Pr! 3 < 1,0 fiir 10? < Re < 10°. Es 
handelt sich dabei um grobe Anndiherung wegen der oben erwihnten von Re abhingigen Faktoren 
und wegen folgender sekundirer Faktoren: Eine untere Grenze fiir die Ubertragung infolge 
Diffusion oder Wiarmeleitung, Konvektion durch Dichteunterschied, sofern nicht Gr < Re?, 
turbulenter und/oder nicht einheitlicher Freistrom, Kinwirkung der Kugelhalterung, Einfluss 
der Kanalwinde, Unsicherheit in den Stoffgréssen der Fliissigkeit, Abhingigkeit der Stoffgréssen 
von der Konzentration oder der Temperatur, Staubansammlung auf der Kugel. 


INTRODUCTION 


Tuts paper the transfer predicted 
from the theory of a laminar boundary layer, 
with experimental results, for local heat and mass 
transfer around the front of a sphere which is 
held in a moving fluid. The theory suggests also 
an approximate but very general expression for 
overall transfer, and experimental deviations are 
shown to be due to the influence of many secondary 
factors. 

The transfer of mass or heat from a 
sphere into a fluid flowing at moderate velocities 
has been treated theoretically by Fréssiine [6] 
and AKsEL’rup [1]. 

Fréss.inc [6] calculated the velocity distribution 
and mass transfer* in the laminar boundary layer 
using power series. By a method outlined in 
Appendix 1 below, he calculated numerically the 


com } sares 


solid 


*For convenience we will refer moinly to mass transfer 
although analogous statements apply to heat transfer. 


coeflicients in the series for the particular case 
of Prandtl numbert 2.532 (corresponding to 
naphthalene evaporating in air). We will express 
the local rate of mass transfer as the dimension- 
less local Nusselt number Nu (@) defined in terms 
of experimental quantities by 

d 2m 


Nu (@) D(C, —C)\I AM), " 


where @ is the polar angle, d the diameter of the 
sphere, D the diffusion coefficient, m the mass of 
the sphere, 4 the area of the sphere, ¢ the time, C, 
is the concentration at saturation and C, that 
Fréssiinc [6] was able to 


in the free stream. 
show that quite generally Nu (@) 0 Re’? and 
that for the particular case of Pr - co (when the 
transfer boundary layer is thin compared to the 
flow boundary layer) Nu(@) « Pr'’*. Conse- 
quently the quantity Nu (@) Re’? Pr’? which 


tFor detinitions see the end of this paper. 
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we shall designate as the transfer number Tr (@) 
approaches a constant value for any given value 
of @. 

Following the method of Frossuinc [6] we 
obtain in the manner shown in Appendix 1 for 


the sphere at Pr 2.532 


Nu (@) Re’? Pr’ 
1-366 (1 - 


Tr (@) 
0-1837 0? 


0-00696 O* + ...) = 1-366 f, (A) (2) 


Because Frésstinc was unable to eliminate the 
Prandt! number from his equations generally 
the coefficients in equation (2) are functions of Pr. 
From the particular solution, Pr — , one finds 


at the front stagnation point (@ 


Tr (0) = Nu (0) Re! * Pr''* 1.472 (3) 


From coetlicients tabulated by Fréssiinc [6] we 
the front 
range 0.5 < Pr < 100 


stagnation 


calculated the transfer at 
for the 
represent the variation with Prandt! number to 


point and can 


0-5 per cent by 
Tr (0) 


1-478 0-158 Pr" for 0.5 Pr © (4) 


Other approrvimate solutions to the transfer 
number at the stagnation point have been derived 
FROSSLING 


K ROUJILINE 


by Fréssiinc [6] and SipunLKrn [25}. 
[6] following PoLmausen [22] and 
[15], used polynomials to express the velocity 
distributions and from his 


and concentration 


solutions we find that, to 1 per cent, 


Tr (0) = 1.53 — 0-190 Pr'*) 0.1 < Pr (5) 


whilst SinpuLKIN using similar expressions but 
with different integral equations gave values which 


we can represent to + 0-5 per cent by 


Tr (0) = 1-50 — 0.182 Pr'’® 06 < Pr <— 10 (6) 


Thus the theoretical transfer numbers at the front 
stagnation point are in quite reasonable agreement 
when calculated by three different methods. Equa- 
tions (4)-(6) indicate that Tr(0) only varies 
about 15 per cent as Pr varies from o down to 0-5, 

AKSEL RUD [1] has also calculated the variation 
in mass transfer around the sphere from the 
front stagnation point to the separation point 


SUTHERLAND 


for Pr 
method. The following equation is derived from 


> o using the approximate polynomial 


AkseL rup’s published graphs 


Tr (@) 1-523 (1 0.1728 @? 


0.0114 0%) 1-523 f, (0) (7) 














Fic. 1. The theoretical distribution of transfer rates 
with polar angle, taking that at the front stagnation 
point as unity. 
theory of Fréssiine at Pr 
f,(@) 


theory of Akset rup at Pr 


Broken curve : 
Equation (2). 
Full curve : x. S(O) 


Equation (7). 


transfer rates calculated 
[6] for 


The distribution of 
by the exact 
Pr = 2-532 — /, (@) in equation (2)—and by the 
approximate Axse rep [1] for 
Pr oo — f,(@) in equation (7) 
Fig. 1. The agreement between the two curves is 
good up to 60°, Beyond this point /, (@) will 
to the absence of terms 


method by Froéss.inc 
method by 
are plotted in 


not be accurate due 
higher than @* which have not been calculated so 
far. Thus /,(@) and f,(@) are independent of 
Re and Pr. Hence Tr (@) is independent of Re 
and varies only slightly with Pr due to the 
variation in Tr (0) discussed above. 

The above theories do not, so far, predict mass 
transfer at the rear of the sphere, or consequently 
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the total mass transfer, because of separation Tr 0-80 (8) 


of the flow with its profound effect on the boundary ; 
. He assumed the transfer to be constant for 


0, << @ < 180 and the separation angle @, to be 
110°, but the experiments presented below show 


layer. We have confined a comparison of experi- 
ment and theory (Section 4b) to the leading 
surface of a sphere. 

Akset’rup [1] calculated the total transfer 
over the sphere by graphical integration and 


that neither of these assumptions is valid. 
However, if we postulate that in the wake a 
obtained velocity and diffusion boundary layer develop 
ince from the rear stagnation point and obey the same 
form of relations as those at the front but differing 
in the values of the constants, then the overall 


2in dGalv. iron pipe transfer number could again be constant exce pt 


cane insofar as the separation point varies and Tr (0) 


varies with Pr. The available experimental data 


_y 


discussed in Section 4, show that T'r does lie 


within fairly narrow limits. 
| 15in. square constant head 
inlet section 
2. EXPERIMENTAL METHOD 


(a) Water tunnel 





Squore grid of 4eéin. aluminiun Uniform flow for velocities of 3-50 em sec7! 
sheet at > in centres 7 
was achieved over a considerable area of the 


> <6 in><—— a 











working section of 3 in. 3 in. of a vertical water 
tunnel by the design shown in Fig. 2(a). The 


25:1 Contraction water was fed from a constant head. the level 


! 
€ 
oO 
N 


being maintained by a centrifugal pump and 
overflow weir. An essential feature is the diffuser, 


3in square transporent shown in Fig. 2 (b), the remaining vorticity after 

working section with . . 
9 passage through this being removed by an 

Groduated screw odopter opening ’ 

~ u 

velocity control : ai x _ J 7 

mechanism section. The design of the contraction section 


aluminium grid at the lower end of the inlet 


(» was supplied by the Aeronautical Research 





3m. squore to 3:n diam 


round section 





4 2ing feed pipe 
3in. g round section 





ane BRwop 


Brass control valve 
and seot 














Gin g flange 
Streamlined support with 


30:n square discharge 
< 9 perforated distributor cone 


receiver Dox 

















=- 30 > 2in. <- 


> 
' 


| 
| 
| 


To flow diffuser 
PESO NIE Ye in. 0.0. porcelain beads 
s 


Contro! valves x 


| Yoin.@ outlet ; Ie O-3in. porcelain insulating 
{ Centrifugal beods in co. |in. layer 
7 recirculating —_|4in. square galvanized iron ‘Perforated zinc screen 


Chain drive Flexible hose connexion pump diffuser box 


























(a) Main section (b) The diffuser 
Fic. 2. The water tunnel 
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Sein. Woll lin LO. rubber tube 
a 


> 

» Pubber gosket 
* 16 TRI 
T A 





= 











-_ Saeerin 











Ye in O.OX7in Brass tube with end piotes 








| 9 











\ 


tein. 1.0. Copper tube soldered into end plotes 


(b) Support cylinder 


0.032 ness. 








(c) 


Fic. 3. 


(a) Unit for inserting sphere into the tunnel 


Laboratory, Dept. of Supply, and was smoothly 
connected into the polymethylmethacrylate work- 
ing section which had 1-5in. diameter screw 
ports at the sides. These ports carried the mounts 
for test objects. A square to round conversion 
section led to the circular outlet where the flow 
was controlled by a blunt-nosed conical valve. 
Water collected in the receiver box was fed to 
the pump by flexible hose to prevent transmission 
of vibration. The tunnel was supported by a 
suitably braced steel frame. 

Most of the construction was in wood, internal 
surfaces being coated with water-resistant resins. 
It would be more desirable—but more costly—to 


Mounting for sphere 


(b) Support evlinder (c) Sphere support 


construct the tunnel in metal which would avoid 
The tunnel required 
regular cleaning, algal growth being controlled 


swelling and distortion. 


by small additions of copper sulphate. The raw 
water supply 
matter to cause difficulties and it was treated 
alum-lime mixture*, the floc being 
allowed to settle for 24 hr when the supernatant 
was found to be sufliciently clear that it did not 
appreciably scatter a light beam. The water 
temperature was controlled by heating or cooling 
coils to 20-0 + 0-3°C, 


contained sufficient suspended 


with an 


*135 galum and 45 g slaked lime in 300 gal water. 
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(b) Spheres and their support 

Each test sphere was moulded on a fine stainless 
wire attached to a heavier support wire (Fig. 3c) 
which fitted into a cylinder (Fig. 3b) carried, as 
a slide fit, by the adaptor (Fig. 3a) which screwed 
into the port. The cylinder acted as a water 
seal and when being removed escape of water 
was prevented by using a clamp closing the rubber 
tube. The adaptor and cylinder were fitted flush 
with the Solution in the time 
(15 sec) required to insert or remove a sphere 


internal wall. 


was negligible. The sphere was centred in the 
tunnel by projecting an image of it and the 
working section on a ground glass screen. 

To reduce any effect which dust may have on 
transfer, the spheres were rotated rapidly for a 
every 10min. Where 
any sign of dust this treatment led to its removal. 


few seconds there was 

The spheres were of benzoic acid, pressed in a 
mild steel mould heated to 110°C. A mark on 
the support wire served to orient the sphere. 
The sphere profiles were measured before and 
after solution either with a co-ordinate travelling 
microscope or from tracings done by hand at 30 
magnifications using a profile projector. The 
spheres before solution were actually prolate 
spheroids the major axis being 0.005 cm longer 
than the minor. The diameter 
0-958 + 0.003 em = and 
1-286 g cm~. 

The rate of solution of the spheres was a 
maximum at the front, a minimum near the point 
of separation and rose to a smaller maximum at 
From the resulting change in shape 


mean was 


the mean density 


the rear. 
measured as changes in radius, the local and 
overall solution rates were calculated as indicated 
in Sections 3(b) and 3(c) below. The spheres 
and their support wire were also weighed before 
and after solution and the overall solution rate 
calculated as shown in Section 3(b) below. The 


experiments lasted between one and four hours, 
in which time the average diameter was reduced 
to about 90 per cent of the initial value. 


(c) Velocity distribution 

The velocity distribution in the working section 
was determined from photographs of the move- 
ment of aluminium flakes. A 250 W a.c. high 


pressure mercury are emitting 100 flashes sec 
tall 
beaker filled with water acting as a cylindrical 


illuminated the working section through a 


lens. The traces were then photographed using 
a 6in. lens. The magnification was determined 


using a known object in the image plane. 


3. EXPERIMENTAL RESULTS 


(a) Uniformity of flow in the water tunnel 


The flow in at least the central 4 cm tom of 
the working section was uniform to within the 
accuracy of measurement of velocity, viz. + 1 
per cent at 3cmsec~! and + 0-1 per cent at 
54cm sec~!, The turbulence, as indicated by the 
maximum observed lateral displacement from the 
vertical path was less than 0-5 per cent in the 


same zone. At 54cm sec 


departure from vertical path could be observed 


up to 3 per cent 


at 1 em from the wall but the internal core of the 
working section was still satisfactory. 


(b) Over-all solution rates 


The sphere changed both its size and shape 
during the several hours in which it was being 
dissolved. The theory presented above applied 
to the initial solution rate. Fig. 4 shows that 
the average solution rate (m, — m,) ¢t during a 
period of solution ¢, varies linearly with the 


square root of the flow velocity, (or Re'’?). As 





»  gmsec” 


Ss 





(me -™)/1 210 








4 6 

US cm? sec’? 

Fig. 4. 

of a benzoic acid sphere (d 

of water flow velocity at 20°C. The straight line is 
(mo m,)/t = 0-439 107-5 u™%. 


The solution rate (averaged over the period 1) 
0-958 em) as a function 
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shown in Appendix 2 the initial solution rate 
can hence be calculated, assuming the sphere 
remains spherical, as 


! 
” 
Zio 


(9) 


= Me — My 


1/2 2 
dt t=O t Me my 


where m, is the initial mass and m, that after a 
time ¢. The observed change in shape was less 
than 0-1 diameters and the linearity of the points 
in Fig. 4 indicates that its effect on the solution 


rate was negligible. The overall Nusselt number 


for the sphere was obtained by substituting the 


initial surface area and the initial solution rate 
derived from equation (9) in equation (1) in 
The diffusion coeflicient of 
benzoic acid in water at 20°C taken as 
0-832 10-5 em? sec™! (Witke ef al. [28]) and 
the saturation solubility as 3-11 g 1-? (KreveLen 
and Krexens [14]). These give Pr — 1210. The 
mean density of the pellets, 1-286 ¢cm-*, was 


integrated form. 


was 


used to calculate the local and integral rates of 
solution from the change in size. 

The integral solution rates are determined by 
substituting equations (10) and (11) in (1) and 
integrating, where, 


Fr Pm & AV, 
dt 0 At 


where 


AV, = (2/2) M0; (75 — 1,) (To sin @; (11) 


is the loss in volume from a zone of the sphere, 


the vr, are the radii after an imterval of time Af 


Fig. 5. The overall solution rate of benzoic acid 
spheres in water in foreed convection. The linear 
relation is given by Nu = 0-582 Rel? Prl/3, 

@ By direct weighing 


A By integration over the surface 


L.. SUTHERLAND 

at positions 0, (7-5 ils. & 3, Bs eo « See 
and A@; — 15. Slight irregularities in the final 
sphere shape caused the calculated solution rates 
to be less accurate than those obtained by direct 
weighing. The variation of Nusselt number with 
Reynolds number obtained by both methods is 
shown in Fig. 5. The results by direct weighing 
were fitted by least mean squares to an equation 
of the form Nu Pr'* = a Re". The best fit is 


Nu Pr * 0-669 Re?! (12) 


with a coetlicient of correlation of 0-984. The best 


fit assuming that n 1/2 is 


Nu Pr''/* = 0.582 Re'/? (13) 


(c) Local solution rates 


The local values of the Nusselt number, Nu (@), 


were obtained from 


Ar (@) (14) 
Pm At 


o-m 
dA MM 


and equation (1) where Ar (@) is the change in 
radius at angle 0. The results are plotted in Fig. 6. 


Pr! * Tr (@) 


' i 
35 45 6580100 “120140160 180 B00 
e° 





Nu (@) Re 





Fic. 6. The local solution rate of benzoic acid spheres 


in water in foreed convection. Separation of 
laminar Reynolds numbers are as shown. 
The broken Nu (0) Pe ¥/3 Reo l/? = os 


(1-00--0-050 07-0-069 04) and is the best fit to the 


layer, 
curve is 


data. 


The best fit of an equation of the form of equation 
(2) to the mean values of Tr (@) at 0°, 30°, 60 
and 90° is 


Tr (@) = 0-94 (1-00 — 0.050 6? — 0.069 @*) (15) 
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(d) Angle of separation 

After a sphere had been dissolving for some 
time, it developed a marked ridge where the 
solution (@ =~ 90°-100°) 
and which was close to the zone of separation of 


rate was a minimum 
the boundary layer. This separation was detected 
by illuminating the 
parallel light and studying the stream lines visible 


system with accurately 


due to differences in refractive index. The separa- 


tion as observed by us from the position of the 
ridge is fitted by the relation 


0 83 — 660 Re? 


(16) 


as shown in Fig. 7 but nearly as good a fit is 


0, — 83 = 191 Re'* (17) 


Expressions (16) and (17) are asymptotic to 


100 





—te!9i Re 

















10° 
3. 7. The angle of separation, 0,, as a function of 
Reynolds number. 
FrOssuine [5]: @ Garner and Grarron [7 
LINTON and SUTHERLAND (this paper); 


TANEDA [26]. 


0, — 83° for large Reynolds numbers ; the value 
of 83° was found by Face [3] at Re = 1-57 « 105, 
Relation (17) is in good agreement with observa- 
tions in the range 15 < Re < 1000, 
GARNER and [7]. by Garner and 
SuCKLING [8] and by Tanepa [26]. Values derived 
from data obtained by Fréssiinc [5] on the 


made by 
GRAFTON 


evaporation of naphthalene are also shown in 
Fig. 7. It is important to note that the overall 
Nusselt numbers are obtained from the summation 
over the surface of the local values and their 
contributions are proportional to the circum- 
ference at each @. Thus the maximum contribu- 
tion is made at the equator and as separation is 
close to this zone it has a considerable influence on 


the overall rate. 


+. D ISCUSSION 
(a) Requirements of Theory 
theories of transfer considered in 


The 


introduction to this paper assume 


the 


(1) potential flow outside the boundary layer, 
uniform far from the sphere. 

(2) the velocity and transfer boundary layers 
are thin compared to the radius of the 

sphere. This ensures negligible contribution 

by molecular diffusion. 

negligible natural convection. 

freedom from turbulence or other inter- 

ference with flow. 


constant fluid properties. 


Table 1. Summary of experiments on the distribution of transfer rates on the sphere 





Tr (0) 


Transfer Experiment 


Mass 36-7: 23 | 
Mass 200-1760 


Heat 1530-4200 


Mass 1490-7580 
Heat t4-15 » 104 


Heat 87-18 = 104 


Theory 
| (degn. (4)) 


Nu (0°-90°) 
Nu (90° 
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Transfer from a sphere into a fluid in laminar flow 


Of the many observations (Tables 1 and 2) 
reported on transfer from a sphere to a flowing 
fluid few satisfy all To 
compare theory with both our own data and 


these assumptions. 


those of others these conditions are examined 
in detail. 


(1) Potential flow (fluid of zero viscosity) 

The theory of Frésstinc can be applied to 
potential flow with any arbitrary velocity distri- 
bution outside the boundary layer. The potential 
flow around a sphere immersed in an infinite 
fluid is assumed in the analysis given in Section 1 
and Appendix 1. Real flow around the front of 
the sphere changes gradually from viscous flow 
at Re - 
Re 


and 


1 to approximately potential flow at 
> 10° (as can be seen by comparing theoretical 

experimental pressure — distributions 
(ScuiicutTine [24] Fig. 1-9). Thus at Re — 20 
the pressure distribution calculated by Kawacuti 
[11] is intermediate between those for viscous 
flow and potential flow. We achieved an initial 
uniform velocity profile in the water tunnel close 
to the throat of the and 
worked at Reynolds numbers of up to several 


contraction section 
thousand to approximate to potential flow. 
Experiments in which the flow was parabolic 
e.g. Garner and Suckiine [8], would not be 
expected to agree accurately with the theory. 


(2) 
For the theory to apply the velocity and transfer 


boundary layers must be thin compared to the 
radius of the sphere. 


Diffusion, and boundary layer thickness 


The latter ensures that 
transfer by diffusion is negligible and requires 
that Re'* Pr'* » 2. This condition often fails 
to apply, particularly to the evaporation of liquid 
drops (see Table 2). Thus the temperature profile 
determined by Ranz and Marsuacet [23] for a 
typical evaporating drop shows that the thickness 
of the temperature boundary layer at the equator 
was 3-6 drop radii. 


(3) 


When a temperature or concentration gradient 


Density convection 


is high, natural convection currents are superim- 
posed on the forced flow altering the rate of 


transfer. Pirer et al. [20], [21] found experi- 


mentally for natural convection alone from flat 
plates that 


Nu = 2 + 0.5 (Gr Pr) (18) 


the 102 < Gr Pr < 108 
the Grashof number. The first term is for diffusion 


in range where Gr is 
and the second term arises from density convection 
and will only be much less than forced convection 
if 
Gr (19) 
Although based on a flat plate, condition (19) 
should also apply to the sphere. 
Table 2 shows that 7'r increases with decreasing 


Re 


importance of the density convection or diffusion. 


and this is consistent with the increasing 


Garner and Suckiinc [8] observed density 
at Re — 200 but 
the effect would be present to the limit given by 
Such effect 
Kramers’ results [13]. We 


on potassium nitrate crystals and spheres of 


convection in horizontal tlow 


equation (19). an is clear also in 


made observations 


sodium hydroxide which showed considerable 
density convection in horizontal flow up to the 
highest Re used by AkseL rub [1] when he studied 
the former material. Our tests are qualitatively 
in agreement with equation (19). 


(4) 
(i) 
stream or arise from the supports for the sphere. 


Turbulence and other interference 
Turbulence. This may be present in the 
Turbulence also arises in the boundary layer of 
10° which is 
MaIse. and 
Suerwoop [18] show that in the range 2000 < Re 

20,000 turbulence in the free stream of intensity 


the sphere, but only above Re 
outside the range considered by us. 


greater than 5 per cent increases mass transfer. 
Some of the data in Table 2 
Tr at Re > 5000 consistent 
turbulence. For Re - 
laminar in effect if the turbulence is less than 


show an increase in 
with increasing 
10* a stream is probably 
5 per cent. Turbulence of 0-5 to 1-5 per cent 
(which we estimate from their range of critical 
Reynolds numbers) may partly account for the 
high values of Tr obtained by Nenakts et al. 
[30] at Re ~ 105, In those experiments of 
Garner and Grarron [7] and Garner and 
SucKLING [8] in which screens were close used to 
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the test sphere excessive turbulence was probably 
present. 

In some of the experiments of Hsu and Sac 
(9), Garner and Grarron [7], Cary [2] and 
XNenakes ef al. [30), the sphere supports seemed 
excessively large or the sphere either was supported 
at the equator or at the front. All these arrange- 
We determined the 
effect which a single side support had at Re 
between 490 and 1120. 
the support were considerably increased and the 
that 


were shed periodically, 


ments lead to turbulence. 


The transfer rates behind 


wake on side oscillated because vortices 


The wake on the Oppo- 
site side of the support was steady but even 
midway the influence of the support could be 
observed. 

With the support Fig. the 


refractive index streamers, which could be readily 


shown in 8 (c) 


observed in parallel light, trailed away from the 


separation point on both sides showing no visible 


asymmetry or instability due to the support. No 
asymmetry in the solution of the sphere due to 
the presence of the support was evident. 

(ii) Duet walls. 
sphere diameter must always be sufliciently large 


The ratio of duct diameter to 


to avoid interaction of their boundary layers. 
(iii) Dust. 
rate is decreased particularly at the front stagna- 


We have observed that the solution 


tion point by deposition of dust. Consequently 
taken but it 


seems to have escaped the attention of other 


care was to eliminate this error. 


investigators. Unfortunately there is no simple 
way of assessing its importance in measurements 


other than our own. 


(5) Variation of fluid properties 

The properties of the fluid vary with concentra- 
tion or temperature and hence between the sphere 
the Where 


are appreciable average values are often used. 


and free stream. these variations 
In most of the work considered here errors due 


to this variation were negligible (e.g. this paper). 
(b) Comparison of theory and experiment 
(1) Local transfer rates 


The curves showing the variation of transfer 
rates with angle are shown in Fig. 8, and these 


SUTHERLAND 


The local heat and mass transfer numbers 


around spheres. 


Fra. &. 


Hist 
[b0) @ 


Separation 


A Cary [2]; and Saat 
1M): A Nexakes ef al, 
SUTHERLAND $ (this 


clata) : 


FrOssuine [5]: 
Linton «and 
paper) ; pont 
(present broken lines, separation points of 


other investigators 


are obtained by taking the average transfer, at a 
given angle, over the ranges of Reynolds numbers 
stated in Table 1. 


out, many data do not satisfy the assumptions 


As we have already pointed 


of the theory. Those of Fréssuinc and ourselves 
appear to do so, and the data by Hsu and Sacre 
on evaporation at the 5 per cent turbulence level 
On 
XENAKIS worked so close to the critical Reynolds 


are possibly satisfactory. the other hand, 
number for a sphere that free stream turbulence 
was probably excessive. Likewise Cary’s work 
has been criticized on the instrumental side by 
Koroskin [12] and his ratio of duct to sphere 
diameter was too small. 

There is a wide variation in these experimental 
results particularly in Tr (0), the transfer number 
at the front stagnation point. In general the 
value of Tr(0) lies below the theoretical value 
(compare equations (4), (5) or (6) with Fig. 8 
and Table 1). Part of the difference may be due 
to the scatter in results and for any investigator 
this is of the same order as those shown in our 
Fig. 6 and seems to be due to the difficulty of 
the measurements. 

The angular distributions are more easily 
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compared by considering the ratios Tr (@) Tr (0) 
as in Fig. 9. The theoretical curve derived by 
AkseL kup [1] is also shown and it is seen that 
agreement between all 


there is substantial 


investigators and theory. 


Relative distribution of local heat and mass 
rutes on the 


Fic. 9. 
transfer 
The 


front portion of a sphere. 


curve is viven by Aksrivrup’s theory 
(sec Fig. 1). 
Symbols as in Fig. 8. except 


GARNER and SUCKLING [8S 


(2) Over-all transfer rates 

The over-all rate of transfer has been studied 
by many investigators (Table 2). The ranges of 
The value of Tr 
was calculated from these data where possible, 


Re, Pr, Gr, Nu and d are wide. 


average duct velocities were replaced by sphere 


approach velocities and the majority of results 


can be represented by 
1-0 for 100 < Re < 10° and 


<Pr.- 


05 < Tr- 


0-5 10% = (20) 


No closer correlation than this seems possible 
hecause the ratio of the integral transfer rate of 
the front half to the rear half decreases with Re 
(shown in Table 1) due to relatively increased 
transfer at the rear. This greatly influences the 
transfer. If Tr(0) varies with Pr as given in 
equation (4) then presumably Tr will vary in a 


like Table 1 the 


variation is only 10 per cent. 


fashion; however, as shows 
AKSEL RUD* [1] derived the over-all transfer for 
1 and Pr 1. The 
agreement is good even for Kramers’ data where 
! Re 0 and 200 < Pr 


the particular case of Re 


BS80.7 


5. CONCLUSIONS 


We 
secondary factors influencing the experimental 
the 
reasonably the distribution of the rate of transfer 


may conclude that, despite the many 


results, boundary layer theory describes 
of either heat or mass around the leading half 
of a sphere. However further investigation is 
needed on the factors which cause theory and 
observation to deviate considerably in transfer 
The form of the 


relation for the overall transfer rates connecting 


at the front stagnation point. 


Nu, Re and Pr is suggested by the theory, and 
results lead to a useful set 
Nu Re '/? Pr? 10 for 
10°, fora variety of systems. 


various experimental 
of limits, Le. 0-5 - 
the range 10? Re 


icknowledgements—-We wish to acknowledge helpful 
suggestions on the manuscript by Dr. IL. R. C. 


C.S.LR.O. 


Prarr, 


Chemical Engineering Section, 


*Frepianper [4] has given a similar theoretical 


analysis but it seems an inferior approximation to 


AKseL rUD'’s analysis. 

+#Two recent papers, (Garner F. H. and Keey R. B., 
Chem. Engng. Sci. 1958 9 119 and 208) on mass transfer 
in forced convection at low Reynolds numbers and in 
available 
The 
experimental results are similar to those presented in 
The work particularly 


natural convection respectively, only became 


to us after the present paper had been written. 


Fig. 6 and Fig. 8 of this paper. 
contirms the important effect of density convection on 
both local and overall transfer rates in forced convection. 
The present authors do not agree with the derivation of 
the theory given in the first paper. 
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AprenbDIx 1 
Laminar boundary layer calculations of the local 
transfer rates on the sphere 


The velocity and concentration distributions in’ the 
boundary layer are determined by the modified Navier 
Stokes, continuity, and mass-transfer equations. These 
for a body of revolution (about the axis of flow) are 


respectively (SCHLICHTING [24]) 


U (a) 
a 


- 
U(x)‘ 


d(ur) 


ou 


d (vr) 0 


(1.3) 
oy 


where (see Fig. 10) y is the distance normal to the surfac. 
of a point in the fluid, a is the distance on the spher 
surface of the normal from the front stagnation point’ 
r is the distance from the axis of a point on the sphere 
surface, u,v are the velocity components in the boundary 
layer parallel and normal to the surface respectively and 
To solve the above equations, 
the 


>» ~), must be known and assuming this 


C is the concentration. 
U(x), the 
layer (i.e. for y 
to be the theoretical non-viscous 


velocity distribution outside boundary 


* potential ” flow we 


have for the sphere of radius r, 


10. Definition of co-ordinates on the sphere. 


S) usin (1.4) 
7 r 


Expressing this in series form 


’ 5 
U (a) “ use 


we have “) 


U,, /80r°, ete. (1.6) 


Similarly the surface of the sphere can be given in series 
form 


ro sin (4 To) (1.7) 


so thet r) = Irs 1 Gr®o. rs = 1/120 r4y, ete. (1.8) 


Equations (1.1) and 1.2) can be combined by introducing 
a stream function ¢ (7, y) defined by 
ddr 


r da 


(1.9) 


(1.2). Within the boundary 
layer the velocity distribution and hence ¢ is a function 


which satisfies equation 
of both «and y and Fréssiine [6] showed that in analogy 


to the series (1.5) ¢ can be expressed 


1/2 


| - [Mi Si (mda 


» 
al 
" 


Zug f(y)? + Bugfg(y)a® 4...) (1.10) 


where /, (»). fg (y), ete. are coefficient functions expressing 


1/2 
Zu 
0(™*) 
’ 


Frossiine [6] modified equations (1.1) (1.11) to make 
in effect, 


the dependence on y and 


(1.11) 


them dimensionless and independent of » by, 
substituting 


y* — y Rel? d, 


with » 
nn” 


with n 


e Rel/? 

y* (2u,*)l/2 and v* = v ~ aaa 
U. 

We will 
(1.1*)-(1.11*). 


Equations (1.1*) and (1.2*) can be solved by substituting 


designate the equations thus modified by 


equations (1.4*)-(1.9*) together with the appropriate 
boundary conditions and equating like coeflicients of 2*. 
One obtains a series of ordinary non-linear differential 
equations which can be solved numerically for the coeffi- 


cient functions /,* (*), £5* (y*), hg*(y*) ete. where 


Jg* (n*) = 85° (n*) (rg* 4,* /1y* Ug*) hy* (y*) ete. (1.13, 


The concentration C in the boundary layer will also 
be a function of 2 and y which may be expressed in series 
form 


C7 (1.14*) 


Cy r*s + 


where the coefficients Co, C5. Cy, ete. are functions of 
y* (or »*). When substituting equation (1.14*) in equation 
(1.3*) a series of ordinary differential equations can be 
obtained for a new set of coefficient functions Fo (n*). 
Gy (y*). i, (»*), etc. when one also substitutes (1.11%) and 
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Fo (y*), Cy 
* a 
‘3 — I, (n*) ete. (1.15*) 


2u.* 
3G, (*) 
es * ry* us* 


") 
The equations involve coeflicients f,*, f,*, ete. and were 
solved numerically by Fréssuinc up to terms involving 24. 

Now the local Nusselt number is given theoretically 
in terms of the concentration gradient in the boundary 
laver by 


Nu (@) 


(1.16) 


Substituting for (9C*  )»*) in equation (1.16) from equation 


ll 
dn* ],e—0 dy* 


(1.14*) with appropriate substitution of the coefficient 
functions from equation (1.15*) we obtain 


Nu (0) Re“? = [— F 


iM, ) re? 


'o (2u,*)'/2 — 2 (2/u,%)*/2 u,* 
(G, 3(2 u,*)! 2 u,* F,’24 —...] 


»* (1.17) 


where dashes denote differentiation with respect to »*. 
For the particular value of Pr — 2-532, corresponding 
and »* = 0, 


W212, 


into air, 
OTH0, F 


101 so that 


to naphthalene 
Fréss.ine [6] calculated Fy 
G, O319, and I, 


vaporization 
4 
1-369 0%? 


Nu (0) Re 2 — 1-862 


O-2705 v*4 (1.18) 


Pr! 3 
02 


sides (remembering 


v* 


both 
inserting 


which on dividing 
Pr = 2-582) and 


equation (2) of the text. 


by 


finally becomes 


APPENDIX 2 


Initial solution rate 


The over-all rate of solution for a sphere is 


dm 


kA(C 
dt . 


C,) 


aD 
d 


Rel/? prl/3 


A and d 


m Wwe have 


as is borne out by the experiments. Expressing 


throughout in terms of m and the density, p 


a 
dt 


(2.3) 


. , 12 
k,, <aPr us (Sx UP) (C, -C.) 


Pm 
Integrating (2.3) 


» 


1/2 
(my*/ 
t 


1/2 
m,*/*) me 


and the initial solution rate is 


. 1/2 
Ky, Me 
t=0 


dm 


dt 


J, (@). 


Sy). 


S*,(y*) 


ry (»*), 


So (0) 


Sg(n) 


» S*, (n*). 


Vr, (n*), 


v* 


az 


F in 


NOTATION 
constant 
area of sphere 
concentration in solution 
cic, ec 
diameter of sphere 


~) 

em 
diffusion coeflicient em? sec~! 
functions of 0 defined in equations (2) and 
(7) respectively 

functions of » defined in Appendix 1 
functions of »* 


£*,(7*), hg* (y* : : , 
star Mg ta? defined in Appendix 1 


: functions of »* 

Gy(9*), My (n°) detined in eceiilid 1 

gravitational acceleration em sec? 

constant defined in Appendix 2 

mass of sphere 

initial mass of sphere 

mass of sphere at a time, f 

radius of a sphere at zero time em 

radius of a sphere zone after a time ¢ em 

constants defined in Appendix 1 

@-'r n= 1,38,5,... 

polar radius, defined in Fig. 10 

change in radius at angle @ 

quantity of heat transferred 

time 

temperature 

temperature on sphere surface 

velocity of fluid in the boundary layer 

parallel to the sphere surface — em sec 

u/U, 

constants defined in Appendix 1 

u,d*/U,,n = 1,3,5,... 

velocity of fluid external to the boundary 

layer em sec~! 

velocity of fluid at infinity (free stream) 
em sec 

velocity of fluid in the boundary layer 

em sec~! 


normal to the sphere surface 
v Re/2/U,, 

distance on the sphere, from the front 
stagnation point, defined in Fig. 10 
aad 


em 
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distance into the fluid, normal to 
the sphere surface, defined in Fig.10 em 
4 Rel? d 
2 2m Local mass transfer 
thermal diffusivity of fluid cm? see~! | ; : : — 
c) 


dm Over-all mass transfer 
C.,) dt coetlicient 


‘ dAM coeflicient 


volume coeflicient of thermal expansion a} 


of thaid per 4 


i/2 Over-all Nusselt number 
y (24, ype 


y* (2u,* yi {2 
thermal conductivity of fluid Nu (@) Local Nusselt number 


cal em! seen) °C-1 


, or: . 
polar angle radian Vr Prandt! number 
angular positions of zones on the sphere Re Udy Revnolds number 
angle subtended by a zone of the sphere For heat transfer 


angle of separation degree 
k , , gd* 8(T, — To) ' 

nematic Viscosity cm sec Gr > Grashof number 
density of thuid gem3 ™ 


density of sphere gem4 dy Over-all heat transfer 
a stream function defined by equation _ dl T,) dt coetlicient 
(1-9) Appendix 1 


Asterisk denotes dimensionless quantities h (0) 


1 24 Local heat transfer 
P coetlicient 


r-7F | dA 
and equations ; (f, Te)\ 

Nu = hdd Over-all Nusselt number 

Subscripts Nu (0) = h(@)d/A Local Nusselt number 


° . » » . . 
Unless otherwise specified Pr = v/a Prandtl number 
rt) refers to initial values (/ 0) For heat and mass transfer 


‘ refers to saturated solution Tr Nu Re-V¥2 Pr-Y¥3 — Over-all transfer number 
f refers to values in the free stream (at intinity) Tr(@) Nu (@) Re V2 Pr-V3 Local transfer number 
0 refers to local values at 0 T 
wits . . ts ransfer number at 
rr (0) = Nu (0) Reo V2 Pre V4 
front stagnation point 

Definitions 
Average Nusselt num- 
bers for 0? < @ mm), 
and 9) <@ 180) 


respectively 


For mass transfer 
ad® (» , Nu (0 mr), Nu (90 180°) 


Gr Grashof number 
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Book Reviews 


Grundlagen der Chemischen Technik. (Herausgeber. 
Hi. Mohler) Fuchs) Band IV —R 
Warmeaustauscher. Verlag Sauerlinder. 
Frankfurt (Main). 500s 45 DM. 


und ©. GREGORIG : 


\arau und 


Das Werk behandelt in seinem ersten Teil den Wiarmeaus- 
tausch im allgemeinen sowie die Probleme des Druckver- 
lustes in Wirmeaustauschapparaten. Anschliessend folgt 
ein kiirzeres Teilkapitel Giber Wirmedurchgang und die 
restlichen dire Teile des Buches sind der wirtschaftlichen 
Konstruktion 
Betrieb von Wirmeaustauschern gewidmet. 
Aufbau des Buches ist 
Das Kapitel iiber wirtschaftliche Dimen- 


Dimensionierung, und dem  technischen 


Der generelle im allvemeinen 
recht gliicklich. 
sionicrung von Warmeaustauschern ist sehr zu begriissen, 
wie auch der Teil iber Konstruktion dem Leser wertvolle 
Unterlagen vermittelt. Dass dabei von der wirtschaftlichen 
Dimensionierung von Isolationen bis zur Festigkeitsberech- 
nung von Dehnungseclementen fast simtliche wichtigen 
Berechnungen cingeschlossen sind, steigert die Nutzan- 
wendung der Monographie fir den Praktiker betriichtlich. 

S>» kurz und sich auf das Wesentliche beschrinkend die 
Kapitel der erzwungenen Stromung abgefasst sind, so weit 
ausholend hat \utor mit Kapitel tiber 
Kondensation befasst, welches fast allzureich mit theore- 


sich der dem 
tischen Betrachtungen von Detailproblemen durchsetzt ist. 

Lobenswert sei hervorgehoben, das fiir viele Kapitel 
Berechnungsbeispicle aufgefiihrt sind, die viel zur Hlu- 
stration der Materie beitragen. Thre Zahl diirfte im beson- 
deren im Kapitel iber Warmetibergang ohne Bedenken 
noch etwas erhéht werden. Dass der Autor in seinem 
Buche, das vor allem fiir den Proktiker gedecht ist. eine 
yudem noch 


so vollstindige Literaturiibersicht, die 


kapitelweise geordnet ist, beifiigt, sei mit Genugtuung 
anerkannt. 

Das Buch kann nicht nur dem praktisch titigen Ingenieur 
und Chemiker in der Industrie empfohlen werden, sondern 
es dirfte auch manchem Studierenden cin willkommenes 
Hilfsmittel sein. 


A. GuYER 


Id. A, Stoff- and Warmeiiber- 
tragung in der chemischen Kinetik. Springer-Verlag 
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FRANK-KAMENETZKI: 
Berlin Gottingen Heidelberg 1959. ix 
DM. 


Es HaNpeLT sich um die deutsche Ubersetzung eines 
1047 in russischer Sprache erschienenen Werkes, von dem 
bereits 1955 eine englische Cbersetzung herausgekommen 
Inhalt ist Zusammenwirken der Trans- 
Stoff mit 


Nach ciner kurzen Einfiihrung in die chemi- 


ist. Sem das 


portvorgiinge von Wiirme und chemischen 


Reaktionen, 
Kinetik 


zahlreiche 


sche und die Wirme- und Stoffiibertragung 


werden Anwendunygscebiete der * Diffusions- 


kinetik 


tet, so 


und der“ 
zB. 


Kohleverbrennung., 


chemischen Hydrodynamik ~ betrach- 
heterogene  Katalyse, 
Ziind- 

Nuch der 


Losevorginge. 
Flammenausbreitung, und 


I Oschvorginge sowhe | eriodische Prozesse. 


sehundire Diffusicnsstrom nach Stefan sowie die Therme- 
diffusion werden berticksichtigt. Der Schwerpunkt der 
Darstellung liegt im Ausbau der Theorie, wobei hiutig 
Ahnlichkeitslehre 


werden 


von der Gebrauch gemacht wurde, 


Seweit vorhanden, Versuchsergebnisse = zum 


Vergleich herangezogen, 


Wir heben es hier mit einem Werk zu tun, das ahnlich 
wie die Monographie von Damkohler (1937) cinen Vorstoss 
in eine noch Disziplin der Naturwissenschaften 
darstellt. 
voller Entwichklung befindet und dem man gerade in den 
letzten Der 


Ubersetzer (Pawlowski Leverkiisen) ist von dem russischen 


junge 
ks wird cin Gebiet abgesteckt, das sich noch in 
Jahren erhohte Aufmerksamkeit zuwendet. 


Original sausgegangen und hat einige Stellen sachkundig 


iberarbeitet und Erginzungen hinzugefiigt. Der deutsch- 


. . 
sprechende Leser wird das Erscheinen dieser Cbersetzung 


beyrtissen. Ihr ist weite Veibreitung in allen interes- 


sierten Kreisen zu wiinschen. 


tl. Gaieuns 
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